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Abstract 

We consider shape and topology optimization for fluids which are governed by the 
Navier-Stokes equations. Shapes are modelled with the help of a phase field approach 
and the solid body is relaxed to be a porous medium. The phase field method uses 
a Ginzburg-Landau functional in order to approximate a perimeter penalization. We 
focus on surface functionals and carefully introduce a new modelling variant, show 
existence of minimizers and derive first order necessary conditions. These conditions 
are related to classical shape derivatives by identifying the sharp interface limit with 
the help of formally matched asymptotic expansions. Finally, we present numerical 
computations based on a Cahn-Hilliard type gradient descent which demonstrate that 
the method can be used to solve shape optimization problems for fluids with the help 
of the new approach. 

Key words. Shape optimization, phase-field method, lift, drag, Navier-Stokes equa¬ 
tions. 

AMS subject classification. 49Q10, 49Q12, 35Q35, 35R35. 

1 Introduction 

Shape optimization problems are a very challenging field in mathematical analysis and 
has attracted more and more attention in the last decade. One of the most discussed and 
oldest problems is certainly the task of finding the shape of a body inside a fluid having 
the least resistance. This problem dates back at least to Newton, who proposed this topic 
in a rotationally symmetric setting. Nowadays, there are a lot of important industrial 
applications leading to this kind of questions. Among others we mention in particular the 
problem of optimizing the shape of airplanes, cars and wind turbine blades in order to have 
least resistance or biomechanical applications like bypass constructions. The wide fields of 
applications may be one of the reasons that shape optimization problems in fluids received 
growing attention recently. Nevertheless, those problems turn out to be very challenging 
and so far no overall mathematical concept has been successful in a general sense. 
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One of the main difficulties certainly is that shape optimization problems are often 
not well-posed, i.e., no minimizer exists, compare for instance |20l [231 EH]. There are 
some contributions leading to mathematically well-posed problem formulations, see for 
instance |25j . but the geometric restrictions are difficult to handle numerically. The most 
common approaches used in practice parametrize the boundary of the unknown optimal 
shape by functions, see for instance mm- However, those formulations do not inherit 
a minimizer in general. For numerical simulations typically shape sensitivity analysis is 
used. Here, one uses local boundary variations in order to find a gradient of the cost 
function with respect to the design variable, which is in this case the shape of the body. 
The necessary calculations are carried out without considering the existence or regularity 
of a minimizer. But in the end one obtains a mathematical structure that can be used for 
numerical implementations. 

In |14j . a phase field approach was introduced for minimizing general volume func¬ 
tionals in a Navier-Stokes flow. For this purpose, the porous medium approach proposed 
by Borrvall and Petersson [I] and a Ginzburg-Landau regularization as in the work of 
Bourdin and Chambolle [5| were combined. The latter is a diffuse interface approximation 
of a perimeter regularization. This leads to a model where existence of a minimizer can be 
guaranteed, and at the same time necessary optimality conditions can be derived and used 
for numerical simulations, see m- In particular, this approach replaces the free boundary 
F of the body B by a diffuse interface. Hence, it is a priori not clear how to deal with 
objective functionals that are defined on the free boundary F. 

In this work, we study the following boundary objective functional: 

J^h{x,Vu,p,u) , (1.1) 

where /i is a given function, u denotes the velocity field of the fluid, p denotes the pressure, 
iz is the inner unit normal of the fluid region, i.e., pointing from the body B into the 
complementary fluid region E = B'^. The velocity u and pressure p are assumed to obey 
the stationary Navier-Stokes equations inside the fluid region E, and the no-slip condition 
on F, namely. 


- div a + {u- V)u = f in E, 

(1.2a) 

divu = 0 in B, 

(1.2b) 

u = 0 on F, 

(1.2c) 


where cr := p(^S/u + (Vti)^) - pi denotes the stress tensor of the velocity field u, p > 0 
denotes the viscosity of the fluid, f denotes an external body force, and I denotes the 
identity tensor. 

An important example of h is the hydrodynamic force component acting on F with 
the force direction defined by the unit vector a: 

h{x, Vu,p, u) = a - (cru) = a ■ + (Vu)^) -pl)iy, (1.3) 


and so (1.1) becomes 


J a ■ {(Tu) d'H'^ ^ crizdB'^ 


(1.4) 


If a is parallel to the direction of the flow, then (1.4) represents the drag of the object 
B. If a is perpendicular to the direction of the flow, then (1.4) represents the lift of the 
object. 
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In the work at hand we propose an approach on how to deal with boundary objective 
functionals in the phase field setting. To be precise, we aim to minimize an appropriate 


phase field approximation of the functional (1.1), and also the functional (1.4), which can 


be considered as one of the most important objectives in shape optimization in fluids. 
The fluid is assumed to be an incompressible, Newtonian fluid described by the stationary 


Navier-Stokes equations (1.2). 


For this purpose, we first discuss how we model the integral over the free boundary T 
if it is replaced by a diffuse interface and how the normal v can be defined in this setting, 
see Section]^ Afterwards, we analyze the phase field problem for both 0 and ( |1.4[ ) and 
discuss the existence of a minimizer and optimality conditions, see Section In Section 
we focus on the hydrodynamic force functional (1.4) and the corresponding phase field 
problem is then related to the sharp interface free boundary problem with a perimeter 
regularization by the method of matched formal asymptotic expansions. We find that the 
formal sharp interface limit of the optimality system gives the same results as can be found 
in the shape sensitivity literature. 

We then solve the phase field problem numerically, see Section For this purpose, 
we derive a gradient flow equation for the reduced objective functional and arrive in a 
Cahn-Hilliard type system. After time discretization, this system is treated in every time 
step by a Newton method. We numerically solve shape optimization problems involving 
drag and the lift-to-drag ratio. 


2 Notation and problem formulation 

Let us assume that II c M'^, d e {2,3}, is a fixed domain with Lipschitz boundary. Inside 
this fixed domain H we may have certain parts filled with fluid, denoted by E, and the 
complement B ■= E is some non-permeable medium. In the following we will denote by 
iz the outer unit normal of B, i.e., the inner unit normal of the fluid region. The aim is to 
minimize the functional, given by (1.1), where F := dB n H, subject to the Navier-Stokes 


equations (1.2). We additionally impose a volume constraint on the amount of fluid. For 
this purpose we choose /? e (-1,1) and only use fluid regions E ci} fulfilling the constraint 

We prescribe some inflow or outflow regions on the boundary of II and choose for this 
purpose g e iT 2(511) such that Jg^g ■ I'on = 0. Additionally, we may have some 

body force / e L‘^{Vt) acting on the design domain. Note that throughout this paper we 
denote valued functions and spaces consisting of M'^-valued functions in boldface. 

As already mentioned in the introduction, problems like this are generally not well- 
posed in the sense that the existence of a minimizer can not be guaranteed. Hence, we 
use an additional perimeter regularization. For this purpose, we add a multiple of the 
perimeter of the obstacle to the cost functional (1.1). In order to properly formulate 
the resulting problem we introduce a design function (/? : H ^ {±1}, where {(p = 1} = E 
describes the fluid region and {(^ = -1} = H is its complement. The volume constraint 
reads in this setting as J^pdx = /3|II|. 

The design functions are chosen to be functions of bounded variation, such that the 
fluid region has finite perimeter, i.e., p e BV{Q,{±1}). We shall write P^{E) for the 
perimeter of some set of bounded variation FI c H in H. Besides, if y? is a function of 
bounded variation, its distributional derivative Jdp is a finite Radon measure and we can 
define the total variation by |D(^| (H). For p e HI/(II, {±1}), it holds that 


|D(^|(II) = 2Pn({</^=l}). 


( 2 . 1 ) 
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For a more detailed introduction to the theory of sets of finite perimeter and functions 
of bounded variation we refer to m m- We hence arrive in the following space of 
admissible design functions: 

$0, := € BV{n, {±1}) I ^ dx = /3 |r!|} . (2.2) 

Let 7 > 0 denote the weighting factor for the perimeter regularization. Then, we 
arrive at the following shape optimization problem for the functional 0 with additional 
perimeter regularization: 


J r 1 7 

-h{x, \/u,p, u^) d iDy^l + - |D(/9| (O), (2.3) 

Q. Z Z 

subject to e and (u,p) e H^{E) x L‘^{E) fulfilling 

-pAu + {u ■ V)u + Vp = f in E = {ip = 1}, (2.4a) 

divu = 0 in FI, (2-4b) 

u = g on n dE, (2-4c) 

u = 0 onr = f7n dE. (2-4d) 


Here, we used the relation (2.1) to replace the perimeter of E with | |D(/5| (H). Fur¬ 
thermore, by the polar decomposition 


D(^ = |D(/?| for e H1/(H, {±1}), 


(2.5) 


of the Radon measure Tip into a positive measure |D(^| and a ^-valued function e 
(H, |D(/?|)'^, see for instance [U Corollary 1.29], we replace the product of the normal 
and the Hausdorff measure in (1.4) by |D(/j|. In particular, Vip can be considered as 

a generalised unit normal on dE. 


We remark that the shape optimization problem (2.3) for the hydrodynamic force 
component (1.3) have been studied extensively in the literature. In the work of [2], the 
boundary integral (1.4) is transformed into a volume integral. This is also done in [71 
[25], but in the latter, the compressible Navier-Stokes equations are considered. We also 
mention |2T], which utilises the approach of Borrvall and Petersson [4] and the volume 
integral formulation. The shape derivatives for general volume and boundary objective 
functionals in Navier-Stokes flow have been derived in [26]. Finally, we mention the work 
of |3], which bears the most similarity to our set-up. Under the assumption that the set 
FI = {<y9 = 1} is and that there is a unique, sufficiently regular solution u to (1.2), the 
analysis of |3| obtained, via the speed method, that the shape derivative of 


J{E) = + 

with respect to vector field V is given by (see O Theorem 4, Equation 39]) 0 

dj(f;)[u] = £{v{o),iy)if-a + pd^q-d.,u)dn'^^\ ( 2 . 6 ) 

^We remark that in [5], the normal n is pointing from the fluid domain to the obstacle, i.e., in comparison 
with our set-up, n = -i/. 
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where q is the solution to the adjoint system (see 0 Equation 33.2]): 

-^Aq + {\/u)’^q - {u ■ \/)q + Vvr = 0 in E, (2.7a) 

divqr = OinE, (2.7b) 

qr = a on r, (2-7c) 

qr = 0 on dQ n dE. (2-7d) 

Here, we denote the normal derivative of a scalar a and of a vector (3 as 

di,a ■■= cx ■ V, du(3-={V(3)v. ( 2 . 8 ) 


We note that as u satisfies the no-slip boundary condition (2.4d), u has no tangential 
components on n dE. Thus, we obtain 


Vw = di,u 0 on r = 0 n dE. 


(2.9) 


Using the divergence free condition (2.4b), and the no-slip condition (2.4d), we obtain on 
T: 


0 = divit = tr{Vu) = ^ d^UiVi = d^u ■ v {yu)^v = {duU ■ u)v = 0, 

i=l 


which in turn implies that 


J{E) = J a-{(Tu)d'H‘^ ^ = J a-{^Vu-pl)v 


m 


d-l 


( 2 . 10 ) 


( 2 . 11 ) 


This is similar to the setting of [261 Remark 12] and by following the computations in |26j 
one obtains (2.7) as the adjoint system and the shape derivative of (2.11) for a domain 
in the direction of U is 0 


DJ(£')[U] = J^{V{0),u) {-pdu{di,u) ■ a + di,p{a ■ u) + iidi,q ■ d^u) dH 

- divr u)^ a - pa) d'H ‘^~^, 


( 2 . 12 ) 


where div r denotes the surface divergence. We introduce the surface gradient of / on T by 
Vr/ with components {Dkf)i<k<d, and with this definition we obtain divr'U = 
for a vector field v. Moreover, in components, we have 


dty{d^u)-a= ^ Uidi{ujdjUk)ak. 
i,j,k=l 

Remark 2.1. In \2fk Remark 12], the term pdu{duu) ■ a appearing on the right hand side 
of (2.12) is originally given as T,ijk=i This is related to di,{di,u) a by the 

formula 


d d^Uk 

y, k'iT. — 7—i^jak = d„{d„u) a- V nidiUjdjUkOk, 
dxidxj 


(2.13) 


where v = denotes an extension of u off the boundary T to a neighbourhood U dT 

with lal = 1 near T and v |r= v. 

^We remark that in |26[ Remark 12] the term divr(li(Vu)a) appears instead of divr(/i(Vii)^a), 
which we believe is a typo. 
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By (2.9), we see that djUk = duU^Vj on F, and so 


d d 

Y, i^idihjdjUkOk = XI T^idihjUjduUkak 
i,j,k=l i,j,k=l 


X l’^idi{\i^j\^)di^Ukak = 0. 
i,j,k=l 


Thus, the last term in (2.13) is zero and we have the relation 


L IZjOk = dr^idr^u) ■ a, 

i,j,k=l (J-kiUXj 


(2.14) 


(2.15) 


when It = 0 on F. 


Based on Remark 2.1 if {u,p) are sufficiently regular, then a short computation in¬ 


volving (2.15) shows that on F, 


- ^divr(( Vit)^a) - fj,di,{di,u) ■ a + di,p{a ■ v) + divr(pa) 

d d 

= Z! i^idk{diUj)ukaj+ Vp-a 

i=l i,j,k=l 

= - pAu ■ a + \Jp ■ a = f ■ a + (u- V)rt ■ a = f ■ a, 


where we have used the no-slip condition (2.4d), and hence (2.12) is equivalent to (2.6). 


3 Derivation of the phase field formulation 


The problem derived in the previous section has several drawbacks. First, it is not clear 
if this is well-posed, i.e., if for every ip e there is a solution of the state equations 
(2.4) and if there exists a minimizer {ip,u,p) of the overall problem (2.3)-(2.4). Second, 
optimizing in the space BV(il.) is not very practical. Deriving optimality conditions is 
not easy and it is not clear how to perform numerical simulations on this problem. Hence, 
we now want to approximate the complex shape optimization problem (2.3)-(2.4) by a 
problem that can be treated by well-known approaches. To this end we introduce a diffuse 
interface version of the free boundary problem by using a phase field approach. 


3.1 The state equations in the phase field setting 

In this setting, the design variable : D M is now allowed to have values in M, instead 
of only the two discrete values ±1, and inherits H^{Q) regularity. In addition to the two 
phases {ip = 1} (fluid region E) and {ip = -1} (solid region B), we also have an interfacial 
region which is related to a small parameter e > 0. By [22], we know that the 

Ginzburg-Landau energy 

T£:iF^(D)^M, T£((^) := ^ I dx + ^V’(v^)dx (3.1) 

approximates ip ^ cq |D(/3| (D) = 2cqP^{{p = 1}) in the sense of F-convergence. Here, 

Co := ^ J^ \/2V'(s)ds (3.2) 

and '0 : M ^ M is a potential with two equal minima at ±1, and in this paper we focus on 
an arbitrary double-well potential satisfying the assumption below: 
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Assumption 3.1. Let ip e be a non-negative function such that ip{s) - 0 if 

and only if s € {±1}, and the following growth condition is fulfilled for some constants 
Cl, C 2 , to > 0 and k > 2: 


Cit^ < 'tp{t) < C2t^ Vt > to- 

Additionally, we use the so-called porous medium approach for the state equations, 
see also [ms]. This means that, we relax the non-permeability of the solid region B 
outside the fluid by placing a porous medium of small permeability « 1 outside the 

fluid region E. In the interfacial region {-1 < < 1} we interpolate between the equations 

describing the flow through the porous medium and the stationary Navier-Stokes equations 
by using an interpolation function satisfying the following assumption: 

Assumption 3.2. We assume that e is non-negative, with ae{l) = 0, ae(-l) = 

cxe > 0, and there exist Sq, s;, e M with Sq < -1 and > 1 such that 

ae{s) = ae{Sa) for S < Sa, 
a£{s) = a£{sb) for s > Sb- 

Moreover, we assume that the inverse permeability vanishes as e \ 0, i.e., lim£^.o«£ = °°- 
In particular, we have that 


0 < ae(s) < sup ae(t) < oo VseM, 

i.e., as € L°°(R). The resulting state equations for the phase field problem are then given 
in the strong form by the following system: 


as{ip)u - pAu + (u ■ \/)u + \Jp = / in n, (3.4a) 

divu = 0 in Q, (3.4b) 

u = g on dCl. (3.4c) 


Later we add |a£((^) 1^1^ dx to the objective functional and this ensures that in the 
limit e \ 0, the velocity u vanishes outside the fluid region, and hence the medium can 
really be considered as non-permeable again. 

In the following, we will use the following function spaces: 


■■= {v e I divu = 0} , := {u e H^{n) \ v\dQ = g, divu = O} , 

and for the pressure we use the space L^{n) ;= {p e L^{Ll) \ f^pdx = O}. The function 
space of admissible design functions for the phase field optimization problem will be given 
correspondingly to (2.2) as 


^ad ■= e | ^ (/?dx = ^ |Q||. 


3.2 The cost functional in the phase field setting 


We are now left to transfer the boundary integral in (2.3) to the diffuse interface setting 
where the free boundary T is replaced by an interfacial region. To this end, we apply a 
result of [22] and approximate the perimeter regularization term with 7^£ s{ip). Mean¬ 
while, keeping in mind the polar decomposition (2.5) and the relation (|2.1|), we consider 
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the vector-valued measure with density as an approximation to ud'H'^ ^. Thus, we 


may approximate (2.3) with 


r 1 7 

/ -h{x,\7u,p,\/ip)dx +—Se{ip). 

Jn 2 2co 

Alternatively, we may appeal to the property of equipartition for the Ginzburg-Landau 
energy, i.e., it holds asymptotically that (see for instance, (5.29) in Section[^ or [9l Section 
5.1]): 


X1 lv»>/ 


dx ~ 0 as £ \ 0. 


Hence, together with ( |2.1[ ), and the fact that T-limit of £s{^) is the functional cq |D(/?| (fl), 
defined for functions with values in {±1}, and +oo otherwise, we have loosely speaking 


2co'H'^-^ L r 


Co iDv^l ~ ^~ 

^ S £ 


(3.5) 


where | and are interpreted as measures on Q, by using their values 

as densities. Here, we have identified r = cl{(^ = l}nf2 with its reduced boundary, then it 
holds that | |D(/?| = |Dx{(^=i}| = l T, see for instance [U Theorem 3.59]. 

The generalised unit normal u can be approximated by To rewrite this into 

a more convenient form, which is in particular differentiable with respect to (p, we use 
equipartition of energy and replace by and obtain the approximation 


CQizd'H 




Tp{ip) 


Vy^dx. 


(3.6) 


Hence, we may also approximate (|2.3l) with 


— f V Vu,p,V(p) dx + 

Co Jn ^ ^ 2co 


(3.7) 


when we extend h{x,Vu,p,-) from unit vectors to all of M" such that h is positively one 

homogeneous with respect to its last variable. This allows us to extract the factor 

We note that in the bulk regions {p = ±1}, we have '4^{p) = 0 and hence the functional 


(3.7) is not differentiable with respect to p. Hence, we add a small constant 5^ to V’ in 


order to have 'ip{s) -i- > 0 for all s e M. However, we neglect the addition of this constant 

for the Ginzburg-Landau regularization £e{^) ia the objective functional because adding 
a constant to the cost functional will not change the optimization problem. 

In fact, for the analysis of the phase field problem, it is only important that 8^ > 0. In 
Section where we perform a formal asymptotic analysis, we will require limes,o<je = 0 at 
a superlinear rate (see Remark 5.1). 


3.3 Optimization problem in the phase field setting 

Combining the above ideas, we arrive in the following phase field approximation: 

1 


min Jg 

(ip,u,p) 


{(p,u,p) := \uf + ^ (I dx 

-I- J M{p)h{x,\/u,p,\/(p)dyi, 


(3.8) 
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subject to e and {u,p) e x Lq^^I) fulfilling 


/ ae{p)u ■ V + fj,Vu ■ Vv + (u ■ V) u ■ V - pdivv dx = / f-vdx '^v€Hq{Q). (3.9) 

*■/ 


Notice, that (3.9) is a weak formulation of the state equations (3.4). Moreover, based 


on the discussions in Section 3.2, the function A4(p) can be chosen to be 

M{ip) = J or M{ip) = — 

2 Co 


(3.10) 


The phase field approximation for the shape optimization problem with the hydrody¬ 
namic force (1.3) is obtained from (3.8) by substituting 

h{x,\/u,p,\/p) = \/ip- {p{\/u + (Vti)^) -pl)a. 

I.e., 


mm 
iV' 


ain Js{(p,u,p):= f ]-ae{p)\u\^ + + -'ip{p)\ dx 

,u,p) Jo. 2 2co \2 e / 

+ J^M.{p)Vp ■ (/u(Vit + (Vw)^) - pi )adx, 
subject to e ^ad and {u,p) e ^.(0) x L^{Q) fulfilling (3.9). 


(3.11) 


3.4 Possible modifications 
3.4.1 Double obstacle potential 

We could also use a double obstacle potential if ): ^ 


)u {-i-oo} instead of the double-well 


potential in Assumption 3.1 


i.e., 






+ 00 


if \ip\ > 1. 


(3.12) 


Then, one has to treat the constraint |(p| < 1 a.e. in the necessary optimality system 
either by writing the gradient equation in form of a variational inequality or by includ¬ 
ing additional Lagrange parameters. Numerical simulations could be implemented by a 
Moreau-Yosida relaxation as in [15]. A Moreau-Yosida relaxation also leads to a differ¬ 
entiable double well potential, and here we restrict ourselves to a differentiable potential 
where both settings can then be included in the above mentioned way. 


3.4.2 Inequality constraint for fluid volume 

Another possible modification of the problem setting would be to replace the equality con¬ 
straint J^ifdx = /3|P| by an inequality constraint /^pdx < /3|P|. This would make sense 
in certain settings, if a maximal amount of fluid that can be used during the optimization 
process is prescribed and not the exact volume fraction. This would not change anything 
in the analysis, only that the Lagrange multiplier for this constraint would have a sign 
and an additional complementarity constraint appears in the optimality system. 
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3.4.3 Objective functionals with no dependency on the unit normal 


We may also consider objective functionals with no dependence on the normal, i.e., the 
boundary objective functional (|2.3|) takes the form 


k{x, \/u,p) 


dn 


d^l 


(3.13) 


An example of (3.13) is the best approximation to a target surface pressure distribution 


in the sense of least squares: 


1 2 

k{x,Vu,p) = - \p-pd\ 


where pd denotes the target surface pressure distribution. Then, using (3.5), we deduce 
that the phase field approximation of (|3.13) is given by 


1 r 1 

— / -ip{ip)k{x,\/u,p) 
Co Jn e 


dx. 


If k{-,-,-) satisfies similar assumptions to Assumptions 4.1 and |4.2| (see below), one can 


adapt the proofs of Theorems |4.6| and 4.10 to obtain existence of a minimiser and the 


corresponding first order necessary optimality conditions. 


4 Analysis of the phase field problem 


In this section we want to analyze the phase held problem (3.8)-(3.9) derived in the 


previous section as a diffuse interface approximation of the shape optimization problem of 


minimizing (1.1) for a Navier-Stokes how. For this purpose, we introduce some notation 


for the nonlinearity in the stationary Navier-Stokes equations. We dehne the trilinear 
form 


b : H\n) X H\n) X H\n) ^ M, 

d ^ 

(u ■ V) V ■ w dx = V / UidiVjWj dx. 

We directly obtain the following properties: 

Lemma 4.1. The form b is well-defined and continuous in the space x H^{Q) x 

Moreover we have: 

\b{u,v,w)\ < A:o||Vu|h 2 (t^)||Vu|h 2 (f^)||Vm||i 2 (f 2 ) Vu,u; e iT(}(S4),u e iT^(S4), (4.1) 

with 

ifd=2, 

^ ifd = 3. 

Additionally, the following properties are satisfied: 

b{u,v,v) = 0 Vu 6 iT^(n), divu = 0, v€Hq{Q), 
b{u,v,w) =-b{u,w,v) Vu e iT^(n), divu = 0, v,w€Hq{^}). 


(4.2) 


(4.3) 

(4.4) 
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Proof. The stated continuity and estimate (4.1) can be found in m Lemma IX. 1.1] and 
(4.3)-(4.4) are considered in [T3l Lemma IX.2.1]. □ 

Besides, we have the following important continuity property: 

Lemma 4.2. Let {Un)neN, {'(^n)neN-, {Wn)nen ■= U,V,W € H^{Ll) be such that Un 

u, Vn ^ V and w in H^{Q) where Vn\dQ. = 'I’bo for all n e N. Then 


hxn. h{un.iVn,w) = b{u.iV,w) 'iw^H^{Pl). 

n^oo 

Moreover, one ean show that 

X 9 {u,v) !->• b{u,-,v) e 

is strongly continuous, and thus 

lim b{un,Vn,Wn) = b{u,v,w). 


(4.5) 


(4.6) 


(4.7) 


Proof. We apply the idea of [321 Lemma 72.5] and make in particular use of the compact 
embedding H^{Q) ^ L^{Ll) and the continuous embedding H^{Ll) ^ L^{Q). The strong 
continuity of (4.6) follows from |32[ Lemma 72.5]. In addition, from the boundedness of 
the sequences (rtn)neN, ('i^n)n6N, (■*^n)n6N, and ( |L6l ), we have 

\b{Un,Vn,Wn) - b{u,V,w)\ 

= \b{Un-U,Vn,Wn)\ + \b{u, Vn, Wn - w)\ + \b{u, Vn - V, w)\ 

< \\Un-u\\L3(n) l|V'Un||i2(0)Hn||L6(0) + l|l^llL6(n)l|Vt7n||L2(0) II“ If ||i3(n) 


<C 


<C 


+ \b{u,Vn-V,w)\ . 


*■0 by MT 


□ 


4.1 Existence results 


In this section, we want to analyze the solvability of the state equations (3.9). Afterwards, 


we will show existence of a minimizer for the overall optimization problem (3.8)-(3.9). 


Lemma 4.3. Let Assumption 3.2 hold. Then, for every e L^{Ll) there exists at least 
one pair {u,p) e Hg ^{Q) x Lq{LI) such that the state equations ( |3.4| ) are fulfilled in the 
sense of (3.9). This solution {u,p) fulfils the estimate 


ll'“lljfi(Q) + bllL2(o) ^ C{p,ae,f,g,n), 


(4.8) 


with a constant C = C{p,ae, f ,g,Ll) independent of p. 

Proof. We refer to m Lemma 4], where the existence and uniqueness statements for the 
velocity field u are discussed. We point out, that the restriction to functions p e L^{Q) 
with l^?] < 1 a.e. in Ll used in |14) is only necessary because the function in jlTj is only 
dehned on the interval [-1,1]. But of course, the same arguments apply to our case where 
as is bounded and p e L^{Ll). 
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Now for every e L^(r2) and u e Hg ^{0.) fulfilling 

y' a^{(p)u ■ V + fj,Vu ■ Vv + {u ■'7)u ■ V dx = J f 


vdx Vnei?'Q^(n), 

we find by 1271 Lemma II.2.1.1] a unique p e LQ{fl) such that ( |3.9[ ) together with 

Ml^q) ^ C{n)\\ae{p)u - pAu + {u-V)u- fWn-^n) 

is fulfilled. Combining this with the previous statements we can conclude the lemma. □ 

This motivates the definition of a set-valued solution operator 

Ss{ip) := {{u,p) e Hg^{Q) X Lo(fl) | {u,p) fulfil ([3^} for p e (4.9) 

Remark 4.1. If there is some u e S^{ip) with || Vu||i;, 2 (f 7 ) < where Kq is defined in 
(4.2). Then Se{p) = {{u,p)}. I.e., there is exactly one solution of (3.9) corresponding to 
ip (see for instance m Lemma 11.5] or m Lemma 5]). 

Moreover, we show a certain continuity property of the solution operator: 


Lemma 4.4. Under Assumption 3.2, assume {ipk)km <= L^{Ll) converges strongly to ip € 
L^{Ll) in the L^-norm and {uk,Pk)km *= x L^{Q) are given such that {uk,Pk) ^ 

Ss{ipk) for all A: e N. Then there is a subsequence, which will he denoted by the same, such 
that {uk,Pk)km converges strongly in x L‘^{Ll) to some element {u,p) e Se{ip). 

Proof. Let {ipk)k(iN and (wfc,PA:)fceN be chosen as in the statement. By passing to another 
subsequence, denoted the same, we can without loss of generality assume that ipk ^ T 
almost everywhere. Invoking (4.8), we obtain a uniform bound on {uk,Pk) in H^{Ll) x 
L?‘{Ll) because {uk,Pk) ^ S£{ipk)- And so there is a subsequence, which will be denoted by 
the same, such that Uk converges weakly in H^{Q) and strongly in L‘^{Q) to some limit 
element u e Hg^^{Q) and pk converges weakly in Lfiifil) to some limit element p e Lq{LI). 

We now aim to show that 


Fk ■■ Hljn) 


f' ^ 2 2 

Fk{v)--= +-|Vu| +{uk-V)uk-v- f-vdx, 

L-converges in ^(0) equipped with the weak topology to 


Foo ■■ Hi (It) 




f' ^ 2 2 

fiv) := + 2 +{u-\/)u-v- f -vdx, 


converging 


as k 00 . To see this we first notice that for any sequence {vk)kGf 
weakly in H^(Q) to n e H^ by Fatou’s lemma it holds that 

/ «£((/?) |u|^ dx < liminf / as{ipk)\vkf dx. 

JQ k^oo JQ 

Applying the boundedness and continuity properties of the trilinear form b(-,-,-), see 
Lemma 4.1 and 4.2, we can deduce that Ihiik^oo b{uk,Uk,Vk) = b{u,u,v). As the re¬ 
maining terms of Fk are weakly lower semicontinuous in and independent of ipk, 

we directly obtain 


Foo{v) < liminf Ffc(ufc). 

k-*oo 
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Let V e be chosen. We will show, that the constant sequence defines a 

recovery sequence. For this purpose, we notice that due to the boundedness and continuity 
of cte, we have from Lebesgue’s dominated convergence theorem 


lim 




v'^ dx 


Invoking (4.5) in Lemma 4. 2l we deduce that 



dx. 


(4.10) 


lim b(uii.,Uk,v) = b{u, u, v), 

k-*oo 


and thus, we obtain that limfe-^tx. Fk{v) - Fao{v). This shows that the F-limit of {Fk)km 
in Hg ^{Q) with respect to the weak topology equals F^o- 

Now we notice, that is exactly the unique minimizer of F^ in as it fulfils 

per dehnition the necessary and sufficient first order optimality conditions for the convex 
optimization problem Fk{u). Hence, the weak H^{Q) limit of {uk)k€N^ which 

is li e Hg ^{Q), is the unique solution of min^^jji Foo{u), thus it holds that 


/ a£{ip)u ■ V + fj,\7u ■ \/v + {u ■ \7)u-vdyi = / f-vdx Vn e 


(4.11) 


By 1271 Lemma 11.2.1.1] we can associate to ( |4.11[ ) a unique p e Lq{^) such that ( |3.9[ ) is 
fulfilled, and hence p - p. Altogether we have shown {u,p) e 

To show the strong convergence in i3'^(n)xL^(fI), we note that from the T-convergence 
of {Fk)keN to Foo we obtain additionally that limfc_^oo Fk{uk) = Foo{u). Invoking Lemma 
4.51 below we find 


lim / a£(v7fc) dx = / ae{p)\u\^ dyi. 
k~*oo JQ JU 


In addition, by means of (4.7) from Lemma 4.2 we have 

lim b{uk,Uk,Uk) = b{u,u,u). 

k^oo 

These two results allow us to deduce from the convergence of the minimal functional values 
of {Fk)km that limfc^oo Jq |Vua:|^ dx = |Vw|^ dx. Then, together with Uk ^ uin H^{Q) 

this yields that limfc„^oo 11^^ - ii|| = 0. 

Subtracting the state equations ( |3.9[ ) written for p from the state equations (3.9) 
written for we find from Lemma |4.5| below that 


/ {pk-p)dwvdy= / {a£{pk)uk-a£{ip)u)-v + pS/{uk-u)-Vvdx 
+ b{uk,Uk,v) -b{u,u,v) 

^ \\o:e{^k)uk - a£(7’)w||i2(Q) H'l’Hiqn) + P ll'^fc “ '“llf/qo) 


k-*oo 
->-0 

+ \\b{uk,uk,-) -b{u,u,-)\\H-i(n) 


k-*oo 
-^0 


-^0 

Thus limfc^oo II V(pfc= 0. Using now the pressure estimate, see for instance [271 
Lemma 11.1.5.4], we find 

\\Pk-p\\mn) ^ c||V(pfc -p)||f/-i(o) 0. 

Therefore, we deduce that {pk)keN converges strongly in L^{Q) to p. □ 
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In the previous proof we made use of the following lemma: 


Lemma 4.5. Under Assumption 3.2, assume that for {(pk)km ^ {'U'k)km ^ L^{Q) 

and (/? e u e L‘^{yi), 


lim \\<pk - <f\\L^{Q.) =0) fk-^ ‘f a.e. and lim \\uk - u||i2(f7) = 0. 

K —► OO K —^ OO 

Then it holds that 

lim / as{<fk) \ukf dx = / ae((/?) |itp dx and lim \\as{(pk)uk - cxeiip)u\\L2(n) = 0 . 
k^oo JQ, JU k^oo 

Proof. Using the ideas of m Theorem 5 . 1 ] and m Theorem 1 ] we find that 
/ aE{Pk)\ukf-ae{p)\uf dx\= / ae{(pk) {\ukf - \uf) dx 

+ J^{ae{(pk) - ae{(p)) \uf dx, 

and from e L°°(M) we obtain 

j^ae{ipk)[\uk\ -|^l| ) dx < ||a£||ioo(K)||rifc+ 'ii||L2(o)ll^tfc-It||i2(n)-^ 0. 

Moreover, the uniform bound on yields by Lebesgue’s dominated convergence theorem 

lim / {asi<fk) - Oisip)) \u\^ dx = 0, 

k-*oo JQ 

which combined with the previous step yields the first assertion. 

Using a similar idea we find 

\\ae{<Pk)uk - aeiT)u\\L 2 (n) ^ \\as{(pk){uk - w)||i2(Q) + \\{ae{<Pk) - a£(v^))u||i2(o) 

^ l|a£||L'>»(R)ll^tfc -w||i2(f^) + \\{ae{(pk) - ae{<p))u\\L2^Q-j 0, 

where we applied Lebesgue’s dominated convergence theorem in order to deduce from 
ae e L°°(M) that limk-,oc \\{aei(pk) - OisiT))u\\L2(^Q) = 0. □ 

We make the following assumption regarding h: 

Assumption 4 . 1 . Let h'. Llx x M x ^ M be a Caratheodory function, which fulfils 

1 . h{-, A, s, w) : Q ^ M is measurable for each w e M'’*, s e M, A e and 

2 . h{x, ■, ■, ■) : X M X ^ M is continuous for almost every x e fl. 

Moreover, there exist non-negative functions a e L^(U), 61,62)^3 ^ L°°{Q) such that for 
almost every x € Ll it holds 

\h{x, A, s, w)\ < a{x) + bi{x) \Af + b2{x) |s|^ + bs{x) |m|^ , 
for all w€R^,s€R,A€ R^^^. 

Furthermore, the functional PL ■ x Lf‘{Ut) x R defined as 

PL{u,p,ip) := L M.{(p)h{x, Vu,p, dx, 
satisfy the following properties 
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is bounded from below, and 


(i) n |j/i_,(n)xL2(n)>«i.,d 
(ii) for all ipn ^ ^ in u in H^{Q), pn^ p in L^(Q), it holds that 

n{u,p,p) < limmfn{Un,Pn, y:>n) ■ 


We then obtain the following existence result for (3.8)-(3.9): 


Theorem 4.6. Under Assumptions 3.1, 3.S\ and\4-l\ there exists at least one minimizer 


of the optimal control problem (3.8)-(3.9). 


Proof. We may restrict ourselves to considering p e with p e [sa, Sf,] a.e. in n. In fact, 
we define as in [221 Proof of Proposition 1] for arbitrary p e the truncated functions 
(f := max{sa, min{(^, Sft}} and find Se{(p) < Seif), where £s is defined in (3.1). Moreover, 


by (3.3), we have a£((p) = oi^{ip) and hence also Se{p) = S^{(p). Therefore we obtain 


J£i‘f,u,p) < J^{p,u,p) for all {u,p) e S^ip) = Ssi<p). 


By Assumption 


4.1 


^ If/i (r 2 )xL 2 (o)x$^jj is bounded below by a constant Cq, and so 
■■ ^ad ^ X Lg(n) is bounded from below by a constant Ci. Thus, we can choose 

a minimizing sequence ((^n, «n,Pn)n€N <= ^ad'>^ Hg^^{9) LI{Ul) with {un,Pn) « Se{pn) for 
all n and 


Wm J^{ipn,Un,Pn) = inf J^{tf,u,p)>-oo. 


In particular, from the non-negativity of f) and a^, we see that for p > 0, there exists an 
N such that n> N implies 


^0 + 2^ ll^¥^n||L2('Q^ < Jg ((p,, 


u 


i,Pn) < inf 


Je{P,U,p)+p. 


Thus, {V(pn}neN is bounded uniformly in L^{Vl). Moreover, without loss of generality, 
we may assume that Pn{x) e [sa^Sb] for a.e. x e 12 and every n e N. And so, we deduce 
that {<Pn}n£N is bounded uniformly in H^{Q) nL°°(fl), and we may choose a subsequence 
(ipnk)keN that converges strongly in L^(fl) and pointwise almost everywhere in 12 to some 
limit element <p e 

Using Lemma 4.4 we can deduce that there is a subsequence of (Un^,,Pnk)keN, denoted 
by the same index, such that 


lim \\un^,-u\\Hi(n) = 0, h™ WPn^-ph^n) = 0, (4.12) 

k-*oo k-^oo 

and {u,p) e Se{p). 

From Lemma |4.5| we deduce additionally that 

lim / a^iifuk) \unkf = / ae{(p)\uf dx. (4.13) 

k-*oo JQ Jfl 

As sup;i,gj^ ||'0(<Pnfc) ||L<>°(n) < 0 ° we can use Lebesgue’s dominated convergence theorem 
to deduce limfc_^oo-(/^((p^.^,) dx = dx. Finally, the weak lower semicontinuity of 

H^{Q) 9 (p Hi- |V(p|^ dx yields 

r ^|V(p|^ + -V'((p)dx <liminf T ^ |V(pnJ^ +-V'(<PnJ dx. (4.14) 

JQ Z £ k-*oo Ju Z £ 
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Together with the lower semicontinuity assumption on % from Assumption |4.1[ we 
deduce that 


J^{ip,u,p) J^{ipnk,Unk,Pnk) = inf 


'e KT'nir ) ; yni, ) ~ Jr ^Pi'^iP')) 


and so {(p,u,p) is a minimizer of (3.8)-(3.9). 


□ 


By the same arguments, one can show an analogous existence result for the optimal 
control problem {(|3.9[), (|3.11[)} involving the hydrodynamic force (1.3): 


Theorem 4.7. Under Assumptions \3.i\ and 3.2, there exists at least one minimizer of the 
optimization problem {( |3.9[ ), p.ll[ )} involving the hydrodynamic force (1.3). 

Proof. We will prove the assertion for the choice Ai((p) = and the analogous 

assertion for the choice | follows along the same lines. 

We first show that {Je{ip,u,p) \ p e {u,p) e S'e((/9)} is bounded from below. We 
may restrict ourselves to considering p e ^ad with p e [sa,Sb] a.e. in as in the proof of 
Theorem 14.61 

Now let p e ^ad be arbitrarily chosen with p e [saj^b] for a.e. x e O and choose 
{u,p) € S£{p). From (4.8), we find a constant C 2 > 0 independent of p such that 

IIpIIl2(o) < C 2 . 

By construction, we have 


P^ [Sa,Sb] => ||V'(<y?)llL<>=(n) < C's, 

for some constant 6*3 > 0 independent of p. Then, using Cauchy-Schwarz’s inequality, and 
Young’s inequality we have 

— f • (//(Vw + (Vii)^)-pi) adx 

Cq Ju 

> - {‘2pC2 + C 2 ) > -^11 -C 4 , 

with some constant C^y Q independent of p. The non-negativity of and ip yield that 

Jeip,u,p) > f —\/^^^^Vp-{p{Vu + {Vu)'^) -pl)a+ ^^\Vpf dx 
JQ Cq JCq Z 

- -^4 > -C 4 . 

This shows that {j£(p,u,p) | p e ^ad,{u,p) e is bounded from below. Hence we 

may choose a minimizing sequence {pn,Un,Pn)neN ^ad ^ x Lg(fl) with 


lim j£{pn,Un,Pn) = mf Je{p,U,p) > -OO. 

n-*oo 

As before, we deduce that {v?n}neN is bounded uniformly in H^{Q) nL°°{Q), together 
with Lemma 4.4, we have subsequences (pni^,Unf.,Pnk)keN, that satisfy 


lim \\pn^ - pWl^q) = 0 , lim 

k-*oo k-*oo 


= 0 , lim \\pn^^ -p\\L^(n) = 0 , 

k-*oo 


16 























and {u,p) € S^{ip). 

To deduce that {(p,u,p) is a minimizer of {(|3.9[), (|3.11[)}, we only need to show that 


liminf / ■ (/U (vitn^ + (VWnJ^) -Pn* I) ctdx 

k-*oo JQ 

^ 'y/V’Cv?) + ■ (/i(Vit + (Vti)^) -pl)adx, 


(4.16) 


as the other integrals in ( 3.11[ ) are shown to be weakly lower semicontinuous in the proof 
of Theorem 4.6 We apply now an idea of [22j and define 


Then we see that 


4>{t) := f \/'tp{s) + Jgds, Wrikix) := 4>{(pnk{x)). 
Jsa 


= (p'{(pnk{x))B(pni,{x) = + ^e)^^nk{x). 

By the uniform boundedness of {(pnj,)k€n in H^{Q) nL°°(Q), we find that ))fceN is 

uniformly bounded in L°°(n), and so by the Cauchy-Schwarz inequality, 

(V’(s) + ds j dx 

< sup (V’(s)+4) / - Sal^ dx, 

se[sa,S6] 

\\^Wnjl2rm < SUp{ll:{ipn^)+5e)\\Dtpnjl2r^y 

fceN 

Thus, we deduce that {wrn.)ken is bounded uniformly in and hence there is a 

subsequence, denoted by the same index, that converges weakly in H^{Q) and pointwise 
almost everywhere in fl to some limit element w e Since (j) is continuous and 

lim^^oo ipnt,{x) = ip{x) for almost every x e 0, we know that w = In particular, the 

weak convergence of DtCnj, to Drc implies that 


^/i/j{iprn:) + SeVipnu ^ +4V(^ in L‘^{Q). (4.17) 

Combining ( 4.12[ ) and (4.17) we obtain from the product of weak-strong convergence: 

lim / \/V’(7’nJ + +(V'Unj'^) - I ) « dx 

_ (4.18) 

= \/'0((p) + 4V(p ■ (/i (vit + (Vw)^) -pi) adx. 

Using ( 4.18[ ), ( |4.13 ) and ( 4.14[ ), we deduce that 

Je(p,u,p) < liminf j£((pnfc,Wnfc,PnJ = inf Je{ip,u,p), 

fc^oo <p<i^ad,{u,p)^Ss(‘P) 

and so {ip,u,p) is a minimizer of {( |3.9| ), ( |3.11[ )}. □ 

Remark 4.2. Note that, for the choice A4((p) = the proof of Theorem ^.1 is completed 
once we showed that Jg is bounded from below, which can be shown similarly as in (4.15), 
and (a) in Assumption jf-Sj has been verified. This follows the product of weak-strong 
convergence: 

lim / Vpnfe ■ (p(Vitnfe + (VWnJ^)-Pnfcl)adx 
kyoJn 

= J Vp ■ (p (Vrt + (Vw)^)-pi) adx. 
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4.2 Optimality conditions 


This section is devoted to the derivation of a first order necessary optimality system for the 
optimal control problem (3.8)-(3.9). For this purpose, we first show Frechet differentiability 
of the solution operator. We will only be able to show differentiability at certain points 
where the solution to the state equations is unique. Otherwise we cannot apply the implicit 
function theorem in order to deduce the statement. To be precise, we obtain the following 
result: 


Lemma 4.8. Under Assumption 3.2, let e nL°°(0) be given such that there 

is with ||Vu£||i 2 (Q) < Then there is a neighborhood N of in 

n L°°(n) such that for every ip € N the solution operator consists of exactly one 
pair, and hence we may write ■ N c H^{^) n ->• x This mapping 

is then differentiable at ipe with DSe{iPe){p>) ='■ {u,p) e Hq{Q) x Lq{Q) being the unique 
solution of the linearized state system 


a^{ipf)LpUe + ae{ipe)u - juAu + {u ■ V)it£ + {ue ■ S/)u + Vp = 0 in n, (4.20a) 

divtt = 0 inUl, (4.20b) 

u = Q on dUt. (4.20c) 


Proof. As already mentioned, we want to apply the implicit function theorem to get the 
statements of the lemma. For this purpose, we first note that, by m Lemma IX.4.2], 
there exists a G e i.e., G satisfies 

div G = 0 in n, G |an= 9- 


We define 

F : {H\n)nL^{n)) X H^{n) X Ll{n) ^ iT“^(fI) X Lg(fl), F= (Fi,F2), 


by 

Fi{ip,u,p)v '.= J' as{(p)u ■ V + pX/u -Vv + {u-'\/)u- V -pdivv - / • r; dx 

+ J {u-\/)G ■ V + {G ■ \/)u ■ V + as{ip)G ■ V + p\/G ■ \/v + {G ■ \/)G ■ vd:x., 

F2{(p,u,p) := divri, 

for all V € Hq(Q). 

Hence, F(ip,u - G,p) = 0 if and only if {u,p) e Se{<p). Thus in particular we have 
F{ip^,u^ - G,ps) = 0. Besides, we directly see that the Frechet differential exists 

and is given at {ipc,u^ - G,Pe) as 


D 


{u,p)Fl{Te, - G,Ps){u,p)v = ae{ipe)u ■ V + pXU ■ VV + {u ■ V)Ue ■ V dx. 

+ J {ui; ■ V)u ■ V - pd\w dx, 

^{-a,p)F 2 {Te,Ue - G,pe){u,p) = divU. 


The assumption ||Vrt£||i 2 (f 2 ) < equations (4.1) and (4.3) ensure that 


ae{ips)u ■ V + p\ju ■ S/V + {u ■ V)lt£ ■ V + {Us ■ '\/)u ■ V dx 
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defines a coercive, continuous bilinear form. Hence, we may use the Lax-Milgram theorem 
and standard results for the solvability of the divergence operator, see for instance m 
Lemma II.2.1.1], in order to obtain that Ug - G,Ps) is an isomorphism. 

Next, we want to consider the differentiability of F with respect to its first argument. 
For this purpose, we have to consider '■ L^{Q) ->• L^iyi) as a Nemytskii operator, 
making in particular use of the embedding ^ L®(n). The results in |29l Section 

4.3.3] ensure that ->• L^{Q) defines a Frechet-differentiable Nemytskii operator, 

which follows from the assumption e L°°(M) n We can then conclude directly 

that F is Frechet differentiable with respect to its first argument with 

D^Fi{ip,u- G,p){(f)v = ■ udx, D^F 2 {ip,u- G,p) = 0. 

Additionally, we need that F is Frechet differentiable in a neighborhood of {(p^,Ue,Pe)- 
To show this, we will use m Proposition 4.14], i.e., we show that the partial derivatives are 
continuous in order to conclude that F is Frechet differentiable. Thus let {'Pk,Uk,Pk)km <= 
n L°°(Q)) X Hq{Q) x Lq{Q) be sequences with 


lim \\uk - u||ffi(n) = 0, ^hm \\pk -p\\L^(n) = 0, lim \\ipk - (p\\m(n)nL--(n) = 0- 


fc-J-OO 


k^oo 


As «£ : L 2 (Q) defines a continuous Nemytskii-operator, making additionally use 

of the continuity properties of the trilinear form as stated in Lemma 4.2, we can deduce 
that 


\\^{u,p)F{Pk,Uk,Pk) - 0- 

Moreover, from a'^ e and standard results for Nemytskii operators we find that 
L^iQ) ^ ip ^ ot'e^P) ^ L®(n) is continuous. And thus we also find by direct calculations 
that limfc_oc ||D^F(v9fc,Ufc,pfc) - = 0. Therefore, we 

obtain that F is Frechet differentiable. 

Finally, applying the implicit function theorem, we obtain for ||v^-¥5£||_H'i(r2)nL~(n) « 1 
the existence and uniqueness of a pair {u,p) such that F{ip,u - G,p) = 0, i.e., {u,p) e 
Ss{ip)- This implies the first part of the statement. The second part of the lemma is a 
consequence of the differentiability statement of the implicit function theorem: 

BSe{<Pe) = - {D(u,p)F{P£,Ue - G,Pe)) ^ o D^F{ips,Us - G,Pe), 

which reads in our setting as div u = 0 and 

/ a'^{ips)(pUi; ■ V + ae{ipe)u ■ V + pVu ■ Vv dx 

r ^ (4.21) 

+ j (u ■'\/)u^ ■ V + (u^ ■'\/)u ■ V - pdivv dx = 0 VueiTQ(n). 

□ 


We denote by Dih{x, A, s, w) for i e {1,2,3,4} as the differential of 
n X X M X 9 (x. A, s, w) h{x, A, s, w) 
with respect to the f-th variable, respectively. 
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Assumption 4.2. In addition to Assumption 4-1, assume further that x 
is in for all (A, s, w) e x M x and the partial derivatives 

D 2 /i(x, ■, s, to), D 3 /i(x, A, ■, to), 04 / 1 ( 2 :, A, s, ■) 


h{x, A, s, to) 


exist for all to e s e M, A e and almost all 2 ; e O. Moreover, we assume that 

|Dj/i(x, A, s, to)| < d{x) + bi{x) |A| + b 2 {x) |s| + 63 (x) |to|, for i e {2,3,4}, (4.22) 

for some non-negative oe L^(0), bi,b 2 ,b^ e L°°{M). 

From Assumption |4.2| we see that 

(0^(0))'^"'^ 9 A ^ D 2 /i(-, A,s,to) 6 l2(0), 

L^{H) 9 s D 3 /i(', a, s, to) e L^(r/), 

(L^(r/))'^ 9 to D 4 /i(', a, s, to) 6 L^(r/), 


are well-defined Nemytskii operators for A e (L"^(0))“’'“, s e L^{VL), and to e (L"^(0))“ if 
and only if (4.22) is fulfilled. Moreover, the operator 

X 0 ^( 0 ) X (0^(0))'^ 9 (A,s,to) 


h{-,A,s,w) € L^(n) 


is continuously Frechet differ entiable, 


Next, by Assumption 


3.1 


we have that Dy(\/'0(y) + 5^) is locally Lipschitz 


and thus the Nemytskii operator 

L°°{Q) 9 ip ^ VV’(v?) + 4 € F°°(0) 


is continuously Frechet differentiable. Hence, we find that 

H ■■ X X n L°°(Q) B {u,p,p) J M.{p)h{x,\7u,p,\/p) dx 

is continuously Frechet differentiable and its distributional derivative is given as 
B'H{u,p,p){v,s,r]) = / A4(¥5)(D2/i,D3/i,D4/i) Vr/)dx 

(4.23) 

+ J h{x,\/u,p,\/p)M'{p)r]dx. 


We note that for the choice A4(p) = the second integral on the right hand side of 
(4.23) vanishes as the Frechet derivative of ^ is the zero functional. On the other hand, 

for the choice A4(p) = : the Frechet derivative is given as 


M'{p) 


1 ' 0 '(<^) 

CO 2^j2{'ip{p) + 6^) 


(4.24) 


Before formulating the optimality system we want to discuss the adjoint system. The 
pair of adjoint variables (qejVrg) e Hq{B) x L?‘{yt) is the weak solution of the adjoint 
system, which is given as follows: find {qe-,T^e) ^ Hq{B) x ifill) such that 


(y-e{‘fe){qe - W) -pdlY (V^e + (VQe)^) + {VUe)^qe - {Us ■ V)qe + VvTe 


= -div {M.{p)B 2 h) 

in H, 

(4.25a) 

divQe = -M{p)B2,h + 'de 

in H, 

(4.25b) 

Qe = 0 

on dH, 

(4.25c) 

where B 2 h,B^h are evaluated at {x,Vu^,pi.,Vpe) and 



de--= i M{p)Bsh{x,VUe,Ps,S/pe) dx. 


(4.26) 
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Remark 4.3. The parameter ??£ e M can be interpreted as a Lagrange multiplier for the 
constraint J^pdx - 0. By carrying out the formal Lagrange method as described for 
instance in and appending the mean value condition on the pressure p with some 

Lagrange multiplier to the Lagrangian, one obtains that appears in the adjoint system 
as in (|4.25|). 


The next lemma shows that the system (4.25) is uniquely solvable: 


Lemma 4.9. Let Assumptions 3.1, 3.2, and f.2 hold, and let e H^{Q) nL°°{Q) and 
Ue € such that II be given. Then there exists a unique solution 

pair (QgjTTg) e Hq{LI) x L?‘{Ll) of the adjoint system (4.25). 


Proof. First, we notice that by definition of d^ (4.26), it holds that 


J^M{ip)T>sh{x,\/Us,Pe,\7pe) -dedx = 0 . 


we 


As (fe 6 L°°{Ll), we have A4((/?) e L°°{Ll) for either choices. Thus, by Assumption 4.2 
obtain that Ai{ip)D 3 h e L^{Q). So, from standard results, see for instance |27l Lemma 
II.2.1.1], we deduce the existence of some w e Hq{LI) such that 

divru = -M.{(p)D^h + d^. 


Note that, by the density of C^^{Q) := {u e (C^(n))'^| diva = 0} in Hq^{LI) (see p7l 
Lemma II.2.2.3]), for any v e there exists a sequence such 

that 


ll'i^'' -^llj/i(!^) 0 as n 


oo. 


Thus, for any y e Hq{Q),v e we find that by the commutativity of second 

derivatives, 

f Vl/'(Vu)^dx= lim f Vj/• dx 

Jn n^oo JQ 

E X E yjdjvyan,i Vjdjdivf dx ) (4.27) 

i^j=l i,j=l 

f {y-V)v^-uondn'^-^ - f y-V(divu’^)dx =0. 

Jan Jn 

We define the bilinear form a ■ Hq^^{Q) x Hq^^{Q) ->• (iTg^^(fl))^ by 
i{u,v) := L ae{p£)u ■ V + yX/u ■ (Vn + (Vn)^) + (■ v - {u^ ■V)u-vdx 


= lim 

n^oo 

= lim 

n^oo 


J a£{ip£)u ■ V + yS/u ■ Vu + (S/u^yu ■ v - {u£ ■ \/)u ■ v dx. 


(4.28) 


where we have used (4.27) for u,v € Making use of ||VrAe||i 2 (Q) < (4.1), 

(4.3), and the Poincare inequality, we can establish that a{-,-) is a coercive bilinear form, 
i.e., there exists a constant c{y, |n|) > 0 such that, 

a{u,u) = / a£{p£)\uf‘ + y\\/uf‘ dx+b{u,u£,u)-b{u£,u,u) 

j o V__/ V j 


=0 by Eis 


>0 

> ^11 Vw|| 1 , 2 - AToll Vri|| 1^2 II VW£||i 2 (Q) > c(/i, |n|)||ll|| 
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Meanwhile, the boundedness of the bilinear form in i?g^(n) x can 

be shown using (4.1), the boundedness of Holder’s inequality and the assumption 
II Vw£||i:, 2 (f 2 ) < Thus, by the Lax-Milgram theorem, we obtain a unique q e ^ (Q) 
such that 

a{q,v) = J ae{‘Pe)us-v + M.{ip){T) 2 h-Vv)(hi - a{w^v) VueiTQ,^(H). (4.29) 


We note that the integral terms are well-defined due to Assumption 4.2 and the bounded¬ 
ness of ae- We set q^ := q + w. The existence of vTg e L^(H) follows from standard results, 
see for instance |27l Lemma II.2.2.1]. Thus, {qe,TT^) is the unique weak solution of the 
adjoint system (4.25). □ 


Now we can formulate necessary optimality conditions for our optimal control problem: 

Theorem 4.10. Let {(p£,u^,p^) e nL°°(H)) x x Lq(H) be a minimizer of 

such that II VUe||i 2 (fj) < Then the following optimality system is fulfilled: There exists 
a Lagrange multiplier Ag e M for the integral constraint such that 


[oi'eiTe) Q \Ue 




) + + Ae -H S/Us,Pe, V^^e), C ) 

/ ZCqS /j 


2coe /L2(n) 

+ (M{ipe)'D4h{x,Vu„p„Vps) + ^yTe,yc) =0 yC^H\n)nL°°{n). 

\ 2co /l 2 (Q) 

(4.30) 

Here, {q^T^e) ^ x L^{Q) is the unique weak solution of the adjoint system (4.25). 


Proof. We rewrite the problem (3.8 )-(|3.9[) as a minimizing problem for a reduced objective 


functional defined on an open set in n L°°(H) by making use of Lemma 4.8 


In 


particular, at least in a neighborhood N c nL°°(H) of pe, the solution operator Ss 

is not set-valued, but for every p € N we have Ss{p) = {(it,p)}. Thus we may dehne the 
reduced functional M by 


j^ip) ■.= J^{p,S,{p)). 

Then, p^ is also a local minimizer of j^. Hence, the gradient equation 

Dj'e (7’e)(7>) = 0, yp€ J^pdii =0, 


(4.31) 


would be fulfilled if would be differentiable. 

We will show in the next step that is differentiable at p^ as a mapping from H^{Ll)n 
L°°{Ll) to M. Lemma |4.8| already ensures that the solution operator is differentiable 
from n L°°(H) to H^{Q) x L^{Q). Thus we now look at dependence of on the 

hrst variable. 

For this purpose we find first as in the proof of Lemma 4.8 that : L^(Q) ->• L^{Q) 
is a Frechet differentiable Nemytskii operator, and hence 

a^ip) \uf‘ dx 

is Frechet differentiable for any u e With similar results, i.e. by making use of 

Section 4.3.3], we also find that 

L°° {Ll) ^ p ^ if{p) e L°° {LI), L°°{Ll)^p^ J' fi{p)d'K, 

H\n) Bp^\/p€ L^(H), H^{n) 3 ^ ^ |V 7 >|^ dx 
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are differentiable. Combining these results and the Frechet differentiability of we find 
that -■ N ^ M. is differentiable. Hence we may conclude by the minimizing property of 


(fi; that the gradient equation (4.31) is fulfilled. We then find from (4.31) that 


0 = = Dj^(v?£)((^) + A£^(/?dx (4.32) 


where we defined 


A, :=-|Hr'Dj,(^,)eM. 


(4.33) 


In particular, we interpret Ag e M as a Lagrange multiplier for the integral constraint 

J^ipdx =/ 3 \n\. 

We now want to rewrite (4.32) into a more convenient form by using the adjoint variable 
Qe, which is defined as the solution of ( |4.25[ ). For this purpose we start calculating the 
derivative of j^. We find for every ip e FP(H) the following formula: 

1 , 


Dj'e (¥?£)</? = \Uef + ae{pe)Ue ’Udx 


'y r 1 

+ — / £Vpe-yp+ -'ll^'{Pe)pdx 
2 co Jo £ 


(4.34) 


+ J^M{pe){^ 2 h,B 3 h,Bih) \(x,vu„p„v<p,) iVu,p,\/p)dx 
+ J^h{x,\/Us,Pe,\/pe)M'{ps)pdx. 

where S^{ps) = {(i^eiPe)} and {u,p) := BS^{pe)p is the solution of the linearized state 
equation (4.20). Now we use the adjoint state as a test function in the linearized state 
equation (4.20) and find that 

a'^{pe)pUe ■ qe + ae{pe)u ■ Qe + pVu ■ VQe dx 
+ qe + {Ue-V)u-qe+p{M{pe)B'ih-'&e) dx =0, 


(4.35) 


(4.36) 


where D 3/1 is evaluated at (x, Vit£,Pe, Vy?e). 

Then we use the linearized state u e as a test function in (4.29) and obtain 

ae{pe)qe ’ U + pV qe ' Vu + (VUej^Qe ■ U - {u^ ■ V)qe ' W dx 
= J^ae{pe)u£ ■ U + M{pe) {B2h ■ Vu) dx, 

where D 2/1 is evaluated at {x,'Vu^,Pe,'^Pe)- 

Comparing (4.35) and (4.36) yields the following identity 

a^{pe)pUe ■ qe + ae{pe)Ue ■ U + M{pe) ( 02 ^ ' Vti +PD 3 / 1 ) dx = 0, (4.37) 

where we have used that p e Lq(H), divUg = 0 in H, u = = 0 on dB, and thus 

f pJedx ='&e f Pdx = 0 , 

Jq Jq 


J^{ue ■ V)q£ ■ w + (Ug ■ v)ii • Qg dx = ^ Wg ■ V(qe ■ u) 


dx = 0 . 
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Hence, by using ( |4.37[ ), we can rewrite (4.34) as follows: 


\ 7^ 


7 

2 coe 


ip'{ipe)(pdx 


Ue\ -Ue-Qe 

+ J^M{ipe)'Dih{x,\/Us,Pe,\/^Pe) ' V if A. 

+ J^h{x,\/Ue,Ps,S/pi;)M'{ip£)(pdx. 

Together with ( |4.32[ ), this yields the statement of the theorem. 

The analogous optimality condition for the optimization problem { ( |3.9[ ) , ( |3.11[ ) } 
volving the hydrodynamic force (m is given as follows: 

Theorem 4.11. Let {ip£,u^,p^) e (<harfnL°°(n)) xiT^ xLq(H) be a minimizer of opti¬ 
mization problem {( |3.9[ ), ( [3.11[ )} involving the hydrodynamic force (1.3) with || < 

thus in particular, Se{pe) = {('“ejPe)}- Then the following optimality system is ful¬ 
filled: There exists a Lagrange multiplier e M for the integral constraint such that 


(4.38) 


□ 


m- 




-\Ue\ -Ue-Qe 


j + + As + M'{ipe)^‘L‘e ' {o-eO) , cj 


L2(r2) 




/L 2 (n) 


= 0 \IC^H^{Ll)nL°°{TL), 


(4.39) 


where cTe '■= ^ (Vw^ + (VUg)^)) -Pel, and {qe,'Ke) ^ Hq{LI) x L^{Ll) is the unique weak 
solution of the adjoint system 

ae{pe){qe “ Ue) - pV ■ {VQe + QeV) + {VUe)^qe - {Ue ■ V)qe + VvTe 

=-/i (div (A4((^£)Vy?e) a - V (A4(:^e)V¥?£) a) inkl, (4.40a) 

dlY Qe = M.{pe)^Pe - U - j- M.{(pe)^Pe ■ udx in Ll, (4.40b) 

Qe = 0 on dkl. (4.40c) 

Proof. Note that for the hydrodynamic force 

h{x,VUe,Pe,V<fs) = ■ {p{VUe + (VWj)^) - Pe I ) ■ a, 

and so we compute that 

D 2/1 = p (Vy?£ ® a + a® Vpe): Ti^h = -a ■ Vpe-i 

B^h = {p{VUe + {VUe)^) -Pel)a. 


As a is a constant vector, (4.22) in Assumption |4.2| is satisfied and the statements follow 
from the application of Theorem |4.10 □ 


Remark 4.4. After using integration by parts, we find that we can rewrite the gradient 
equation (4.39) for the hydrodynamic force formally in the strong form as 

|2 


7 

2 co 


+K + a'eiTs) Q \Ue 


Ue-Qe 


j - M.{tpe) div {(Tea) = 0 m n, 

(4.41) 


with the boundary condition 


1 

—eVipe ■ ^dn ■^M{pe)vdn ■ {(^eo) = 0 on dPL. 
2 co 


(4.42) 
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Moreover, with sujficiently smooth solutions, we can make use of the state equation 


(3.4a) to rewrite (4.41) as: 

(sAipe - + Ae + a'^{pe) Q 


2 Co 


-Ue 


Qe 


(4.43) 


+ Mipe) if - as{Pe)Ue - {u^ ' V)u^) ■a = 0. 


Remark 4.5. We note that the above analysis of (3.8)-(3.9) can be modified to include 
a Dirichlet condition for the design function on dQ, for instance = 1 on dQ. This 
amounts to changing the space of admissible design functions to 

^ad = ^ I J' ip d-K. = /3 \Q,\ and p = 1 on . 

Then, in the optimality conditions ( |4.30[ ) and (4.39), and also in (4.31) and ( 4.32 ), we 
use test functions C e n and p e Hq(Q). Moreover, from Remark f.f, the 

strong form of the resulting gradient equation (4.39) remains as (4.41) (or (4.43)), but 
now with the boundary condition 


Pe = ^ on dTl. 


5 Sharp interface asymptotics for the hydrodynamic force 


In Section]^ we introduced the diffuse interface problem ( |3.8[ )-( [3(^ as an approximation 
of the shape optimization problem (2.3)-(2.4) for a general functional h. In Section]^ 


the existence of a minimizer {p^,u^,p^) to (3.8)-(3.9) for every fixed e > 0 is guaranteed 


by Theorem 4.6, and the first order necessary optimality condition is given in Theorem 
The analogous results for the hydrodynamic force problem {(|3.9[), p.ll[)} are also 


4.10 


presented in Theorem 4.7 and Theorem |4.11 


In this section, we focus only the hydrodynamic force problem {( |3.9[ ), ( |3.11[ )} and carry 
out a sharp interface limit of the system {(3.4), (4.40), ( 4.43[ )} by the method of formally 
matched asymptotic expansions. We hereby recover the optimality conditions expected 
by classical shape sensitivity analysis presented in Section in the limit e \ 0. For an 
introduction and more detailed discussion of the techniques and basic assumptions used 
in the method of formally matched asymptotic analysis we refer for instance to HUES]. 

In the asymptotic analysis, we assume there are sufficient smooth solutions to the 


system {(3.4), (4.40), (4.41)}, and hence we consider (4.43) instead of (4.41) in the sequel 


as the analysis is comparatively easier. 

Assumption 5.1. We assume that for small e, the domain hi can be divided into two 
open subdomains 11'''(e), separated by an interface r(e). Furthermore, we assume that 


there is a family {pe,Ue,Pe,q£,T^£,\e,'&e)e>o of solutions to {(3.4), (4.40), (4.43)}, which 
are sufficiently smooth and have an asymptotic expansion in e in the bulk regions away 


from r(e) (the outer expansion, see Section 5.1), and another expansion in the interfacial 


region (inner expansions, see Section 5.2), see also UMIW for a detailed formulation. 


For the remainder of this section, we will make use of the following assumptions ex¬ 
tensively: 

Assumption 5.2. The correction constant 5^ and the interpolation function fulfill 

de = , k> 1, (y.£{t) = -d{t), 
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(5.1) 


where a e n L°°(M) satisfies the following properties: 

d(-l) > 0, a(l) = a^(l) = 0, a{t) ^ 0 for i ^ 1. 

Moreover, we assume that the potential fi e C'^(M) satisfies: 

= fi\±l) = 0. (5.2) 

For the terms involving the square root, we make use of the following expansion for 
a = ao + eai + 6^02 + ..., which holds due to Taylor’s theorem: 

\/a + Ss = y/ao + eoi + ... + e^{ak + 1) + ... 

+ ...+£: (a^ + 1) + ... j + .... 



5.1 Outer expansions 

We assume that for e {(ps,u^,p^,Xs,'&^,q^,7r^}, the following outer expansions hold: 


Ve = vo + evi + .. 


Applying Taylor’s theorem and (5.3), for the choice A4{(ps) = we obtain 

following outer expansion 


M{ipe) = M{ipQ+£ipi + ...) 
1 


\/2cq 

1 

\/2co 


[\J V’(¥’o) + ^/’'(v?o)(e<y?i +/...) + ... + + ... )^ + ... j 


(5.4) 


2\/V’(¥^o) 


+ 0{e ) =: A4o(v^o) + £Mi{ipo)tpi + h.o.t.. 


We remark that, for the classical smooth double-well potential V’(¥^) = |(1 one can 

compute that 


lim^ 


fi'{s) 
2, hm —= 
\Jfi{s) 


and so A4i(±l) is well-dehned for the smooth double-well potential. 
We denote {■)q t o be t he order fi outer expansions of equation (■). 
To leading order (‘iAa)~Q gives 

a{ipQ)uQ = 0 . 

= 0. Similar!] 

«(</^o)qo = a{ipQ)uQ. 


By (5.1), if ipQ ^ 1, we then obtain uq = 0. Similarly, to leading order (4.400)^^ gives 


(5.5) 

(5.6) 


Thus, if ifo +1, then go = up = 0. 


Meanwhile, (3.46)^, (3.4c)q, and (4.40c)q give 


divrto = 0 in n, 
uq = g, Qo = 0 on dQ. 
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To order -1, {4A3)q gives 




(5.7) 


If (fo 1, then from (5.5), (5.6), and (5.1), we have that 


-'ijj'iipo) = 0. (5.8) 

Hence, (/pq must be a piecewise constant function that takes values equal to the roots of 
ijj'i-). The stable solutions to (5.8) are ipo = ±1. In particular, we can define the fluid 
region and the solid region by 

E ■■= {x € ^}\ ipo{x) = 1}, B ■■= {x €Q \ ipo{x) = -1}, 


respectively. Moreover, from ( |5.5[ ) and (5.6) we have 

uq = Qo = 0 in B. 


(5.9) 


Furthermore, as (pQ = ±1, we have Vpo = 0 in Fi and H, and so, from the definition ( 4.26[ ) 
that = 0. From (4.406 )q we have 


divQo = 0 in Fi u H. 


The next order (3.4a )q gives 

a'{po)piUo + a{Lpo)ui - pAuo + (uq • V)wo + Vpo = /• 


By (5.1), for (po = 1, we obtain 

-fiAuQ + {uq ■ V)uq + Vpo = / in Fi. 


(5.10) 

(5.11) 

(5.12) 


Similarly, (4.40a )q gives 

a'{po)^i{(lo - Uq) + a{po){qi - ui) - ^udiv (Vqo + (Vqo)^) 
+ (Vuo)^qo - (uq ■ V)qo + Vtto = 0. 

For (y9o = 1, we obtain 

-fiAqo + (Vwo)^Qo - {uo ■ V)qo + Vvtq = 0 in .F, 


(5.13) 


where we have used (5.10) to simplify the divergence term. 


5.2 Inner expansions and matching conditions 

Now we consider the interfacial region, i.e. near some free boundary F = dE n dB which 
is assumed to be the limiting hypersurface of the zero level sets of p^- For studying the 
limiting behaviour in these parts of Q we introduce new coordinates. For this purpose we 
introduce the signed distance function d{x) to F and set z = ^ as the rescaled distance 
variable. Here we use the sign convention d{x) > 0 ii x € E. 

Let 7 (s) denote a parametrization of F by arc-length s, and let u denote the outward 
unit normal of F. Then, in a tubular neighbourhood of F, for sufficiently smooth function 
v{x), we have 


v{x) = v{'y{s) + ezu{'y{s))) =■ V{s,z). 
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In this new (s, z)-coordinate system, the following change of variables apply, see m-- 

Vx'W =+ h.o.t., 
e 

where Vr/ denotes the surface gradient of / on T with components {D.kf)i<k<d and h.o.t. 
denotes higher order terms with respect to e. Moreover, if n is a vector-valued function, 
then we obtain 

diVajU = ^dzV ■ V + divpV + h.o.t.. 

e 

In particular, using the fact that the normal v is independent of z, we have 

An = divxCVxw) = \dzzV + - divr(<9^1^^') + h.o.t., 

e V_^_- 

=-KdzV 

where k = -divri^ is the mean curvature. 

We denote the variables ipe, u^, p^, in the new coordinate system by <!>£, Ue, Pe, 

Qe, n^- We further assume that they have the following inner expansions: 


I4(s, z) = 1 / 0 ( 5 , z) + £Vi{s, z) + 
for We then obtain, 

= >lo(^o)+£-Adi(^o)^i + h.o.t.. 


where A^o, ■M.i are as defined in (5.4) if we consider A4(y?) = ■ 

We remark that, for a sufficiently smooth function / independent of e. 


f{x) = fijis) + £zu{s)) = fijis)) + £z\/f{-f{s)) ■ u + h.o.t. 
=■ Fo{s) + £Fi{s, z) + h.o.t., 

for X in a neighbourhood of F. As a consequence, we see that 


dzFo = 0 . 


(5.14) 


As the Lagrange multipliers Ag and 19^ are constant, we assume that the inner expansions 
are the same as the outer expansions. In particular, the leading order expansions of the 
Lagrange multipliers do not depend on z. 

The assumption that the zero level set of converge to F implies that 

4>o(0) = 0. (5.15) 


In order to match the inner expansions valid in the interfacial region to the outer expan¬ 
sions of Section 5.1 we employ the matching conditions, (for the derivation we refer to |16L 
Appendix D]): 


lim Vb(5,z)=no, (5.16) 

Z->±c>o 

lim dzVo{s,z) = 0, (5-17) 

Z->±oo 

lim dzVi{s,z) = Vvn ■ v, (5.18) 

Z->±oo 

lim dzzV 2 {s,z) = {{v'-\7){u-\7)uQ) = diz{dizUQ), (5.19) 

Z^±oo 
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where r;Q ;= lim^^o ^o(p ± for p e F. Then (5.18) and (5.19) for vector-valued functions 
read as 


lim dzVi{s,z) = duVQ, lim d^zViis, z) = ■ V (Oi^Vq) = di^{di^VQ). 

2-»-±00 Z-*-±00 

As div Us = 0, we can rewrite 


Aus = div (Vtle + (VUs)'^). 


For a tensor A, let S{A) = ^{A + A^). Then we can compute 

Aus = ^dz{£{dzUs 0 u)u) + -dz{£{S/rUs)iy) + - divr{£{dzUs 0 iv) + ... 

e e 

= \dzzUe + \dz{dzUs ■ v)v + -dz{£{VvUs)v) + - diYr{£{dzUs ®u) + .... 


We note that the same ex pans ion holds for the divergence term in (4.40a). 


Similarly as in Section 5.1, we will denote (•)j to be the order /? inner expansions of 
equation (•). 


5.2.1 Inner expansions of the state equations 

To order -1, (3.46)J^ gives 

dzUo ■ V = dz{Uo ■u) = 0, 


while to leading order (3.4o)r^ gives 


-ndz{dzUo + {dzUa ■ v)u) = -[idzzUo = 0, 


(5.20) 


(5.21) 


matching condition (5.17) leads to 


where we have used (5.20). Integrating with respect to z from -oo to z and applying the 


dzUois,z) = 0, (5.22) 

and so Uq is independent of z. Integrating once more with respect to z from -oo to z and 


by the matching condition (5.16), we hence find that 

Uq{s,z) E Uq = 0, 


where we made in particular use of (5.9). This implies 


= no = 0. 


To first order (3.46)^ gives 


dzUi ■ V + divri/o = dzU\ • n = 0, 


(5.23) 


(5.24) 


(5.25) 


where we have used (5.23). Using (5.23) and (5.25), to first order (3.4a)j^ gives 


-^idzzUi + dzPou = 0 . 


(5.26) 
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5.2.2 Phase field equation to leading order 

To leading order (4.43)J^ gives 


- V^'(^o)) + d'(^o)(i \Uo? - Uo ■ Qo) - M{^o)ai^o)Uo ■ a = 0 (5.27) 

2co 


Using (5.23), the above simplifies to 

d,,4>o-V''(^o) = 0. 


(5.28) 


Along with the matching conditions (5.16) for 4>o: 

4>o(s,z = ±oo) = ±1, 

we can choose 4>o to be independent of s and as the unique monotone solution to ( 5.28| ) 
satisfying 4>o(z = 0) = 0 (recall (5.15)). Moreover, taking the product of (5.28) with <i>Q( 2 ;) 
and integrating with respect to z from -oo to z leads to the so-called equipartition of 
energy after matching: 


^ |^o(^)f = '^{^o{z)) for |z| < oo. 


(5.29) 


Moreover, a short calculation using (5.29), the monotonicity of 4>o, and a change of vari¬ 
ables s $ 0 ( 2 ) shows that 

CO = ^ ^1 \/2V’(s) ds = ^ ^ \/2V’(^o(z))4>o(z) dz = ^ ^ |^o(^)P dz . (5.30) 


5.2.3 Inner expansions of the adjoint equation 

Before we analyse the adjoint equation, we first compute: 

diY{M{(pe)Vipe) = + divr ^A4(^£) -h 

+ h.o.t.. 


and for any 1 < j < d, 


{\/{M{(pe)\/iPe)a)j = Y,di{M{ipe)djipe)ai 


2=1 


= ^ ai + Di ^A4(4>£) j j a* 


-H h.o.t., 


so that 

S/{M{(pe)V(pe)a = a)udz{M{^e)dz^e) 

+ ^ (( 1 / • a)d^(A4(4>£)Vr^>£) + Vr(A^(^£)d^4>£Zz)a) 
+ VrCVr^heja-H h.o.t.. 


To leading order (|4.406l)j^ gives 


d^Qo • = A4o(^o)^o(*^ • “)> 


(5.32) 


(5.33) 
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while to leading order {AA^a)j gives 


- ^J-^z{^^QQ ■ v)u = -^a^(A1o(^o)^o)((*^ • a)u + a), 


(5.34) 


where we have used (5.31), ( |5.32[ ) and that u is independent of z to simplify the right 
hand side of (4.40a 


Integrating (5.34) with respect to z from -oo to z and using the matching condition 
dzQo + {dzQo ■ I')^’ = A4o(^o)^o((i^ • a)v + a), 


(5.17) leads to 


(5.35) 


and upon adding the product of (5.33) with u leads to 

dzQois, z) = A4o(^o)^o®- 

Integrating (5.35) with respect to z from -oo to z, using the matching condition ( 5.16[ ) 
and Qo = 0 (see ( |5.9[ )), lead to 


Qo{s, z) = A4o(^o(^))^o(2) dz j 


a. 


(5.36) 


In particular, the right hand side is independent of s, and so we can deduce that Qo is 

(5.37) 


also independent of s. Using the matching condition (5.16), we hence have 


9o 


A4o(^o(^))^o(2)dz j 


a. 


For the choice A4((/?) = |, we see that 

J^Mo{^o{z))^q{z) dz = ^ ^>o(^) dz = 1, 

while for the choice ■M((p) = we see that by (5.4), (5.29), and (5.30), 


(5.38) 


A4o($o(2))^o(2) dz = — £ -£=a/V'($o(^))^o(2 ) dz = — £ J |$o(z)f dz = 1. 

Cq -J M Cq M .Z 


Thus, in both cases, we obtain 


Qo = “• 


To the next order, we obtain from (4.406)) 


dzQi ■ iz = MQ{^o)d^^i{u ■ a) + A4i($o)^i^o(*z • a), 


(5.39) 

(5.40) 


where we used that Qq and 4>o are fu nctions of 2 ; only, and 'i?o 
Meanwhile, from (|5.31l) and (|5.32l), (4.40a)7^ gives 


0 from the outer expansions. 


a{^Q)Qo- ^idzzQi- lJidz{dzQi ■ iz)iz - 2/r(92(T(VrQo))iz - 2^divr(T(Qo 0 u)) 

= -/i(a + (iz-a)i/)9^(A4o(^o)92 4>i + A4i(4>o)4>i$o) (5.41) 

- /idivr(A4o(4>o)^oiz)a - /rVr(A^o(^o)^oi^)a- 


Moreover, we can simplify, thanks to fact that <I>o and Qq only depend on z: 


2divr(T(Qo 0 u)) = Vr(Qo)iz + (divriz)Qo + (Vriz)Qo + (divrQo)iz 

= ~hQq + (Vriz)Qg, 

divr(A4o(^o)^o*z) = -A4o(^o)^o«^, 

Vr(A4o(‘^o)‘l’o*z) = A4o(^o)‘l’o^r*z. 
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Then, using the relation (5.35), we obtain from ( |5.41 ): 

a(^o)<3o - - IJ^dzidzQi ■ u)u + ^jlkQ'q - /x(Vr*^)Qo 

= - //(a + {v ■ a)v)d:,{Mo{^Q)dz^i + A 4 i( 4 'o) 4 >i 4 >o) + - ^{Vvi')Q'o, 

and thus, upon cancelling the common terms, we have 

a(^o)Qo - fJ^dzzQi - fJ-dz{dzQi ■ i')i' 

= - /i(a + {u ■ a)v)dz{Mo{^Q)dz^i + A4i(4'o)4'i4'o). 

5.2.4 Phase field equation to first order 

Using (5.23), we obtain from (4.43)5 to first order: 

^ [-dzz^l + + V^"($o)$l) + Ao 

2co 

- a{^o)Ui ■ Qo + A4o(4’o) (To - a{^o)Ui) ■ a = 0. 


(5.42) 


(5.43) 


(5.44) 


(5.45) 


■\Z=+oo 
\z=-oo • 


Making use of (5.35), after taking the product of (5.44) with <hQ we have 

+ 4>i(V''(^o))' + K + Ao$o 

ZCq ' / 

- (d(4>o))'f7i • Qo + Qo • (To - a($o)C/i) = 0. 

We note that by integrating by parts: 

- /'(d($o))'?Ti-Qodz = /'d(4>o)(f7i-Q[, + d,f7i-Qo)dz-[d($o)t/i-Qo]) 

We use that q;(1) = 0, Qo{z = -oo) = Qq = 0 to deduce that the jump term is zero. Hence, 

- f (d(4>o))'C/i-Qodz = [ a{<^o){Ui-Q’o + dzUi-Qo)dz. (5.46) 

So, from integrating ( |5.45[ ) over M and using (5.46) we obtain 

f i-dzz^i^'o + 4>i(V''(^o))' + K l^oH + Ao4>o dz 
Jr zco ' ' ' / 

+ f d(4'o)<9^f7i • Qo + Qo • Todz = 0. 

-J M 


(5.47) 


Considering the first line, we find that, after integrating by parts and applying match¬ 
ing (5.16)-( 5T7l ) for <ho. 


f 7^ i-dzz^i^'o + 4>i(^/^'(4>o))' + K|$of) + Ao4>odz 

JR 2co ' ' ' / 


^ f d,<hi(4>''-V^'(ci>o)) + 7^[V>'(d>o)<hi-ck'd,$i]:: 

2 co Jr 2 co 

f l^of dz+Ao f 4>odz = K 7 + 2 A 0 , 

2 cn Jm. Jr 


(5.48) 


=2co 


where we made use of (5.28), the relation (5.30), and that k is independent of z. Thus it 
remains to identify 


f Fo-dzQo + d{^o)dzUi-Qodz. 

•J M 


(5.49) 
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To this end, we take the scalar product of (5.43) with dzUi and use (5.25) to obtain 


(5.50) 


a($o)<3o • dzUi - ndzzQi ■ dzUi 
= - ^l^zUl ■ adz{Mo{^o)dz^i + A4i($o)$i^o). 

Integrating ( |5.50 ) over M with respect to z, and applying integration by parts leads to 

f a{‘^o)QQ-dzUidz 

= f dzzQi-dzUi-dzUi-adz{Mo{^o)dz^i+Mi{<^o)^i^'o)dz 

Jm 

= IX [dzQi • dzUi - dzUi • a{Mo{^o)dz^i + 

- IX r dzzUi • (9,Qi -a(A4o(^o)5.$i +A^i(^o)$i^o))dz. 

•J M 


(5.51) 


Using ( |5.2[ ), the matching conditions ( |5.16[ ), ( |5.17[ ), ( |5.18[ ) for 4>o, and ( |5.18[ ) for Qi 
and LIi, we see that the jump term is 

[dzQi ■ dzUi - dzUi ■ a(7Wo($o)5.4>i + A4i($o)$i^o)];:!lI = [5.90 ■ 5.«o]!, (5.52) 
since, in the case M.{^) = we have A4o(4>o) = \ and A4i(4>o)4>i = 0, while for the case 
Ai((p) = using (5.29) and the matching conditions, we have 

[dzUi ■ a(A4o(4>o)5,$i + Mi(d>o)^i^o)]niI 
1 


y2co 


dzUi ■ a 


\ n/2 ^2V>(4>o), 


(5.53) 


Meanwhile, using (5.40) and (5.26), the integral term is 

f fidzzUi ■ {dzQi - a(Mo(d>o)5.$i + A4i(4>o)^i$o)) dz 
= f -dzPQf ■ {dzQi - a(A4o(4>o)524>i + A4i(4'o)4>i4>o)) dz = 0. 

JE 

Together with ( |5.14 ), i.e., Tq is independent of 2;, we obtain from (5.52), (5.53) that 
f Fo ■ dzQo + a{^o)dzUi -Qodz = /o • [90]- + Ix[5.9o • 

= /o ■ a + liS.Qo ■ 5.rxo, 
as Qq = Uq = 0, and Qq = a from (5.39). Thus, we obtain from (5.47) the following 
solvability condition for 4>i: 

2Ao + K7 + /o ■ a + IxS.Qq ■ S.rtg = 0 on T. 


(5.49) is 


(5.54) 


5.2.5 Sharp interface limit 

In summary, we obtain the following sharp interface limit: 

-fj,Auo + (uq ■ V)wo + Vpo = / in Til, (5.55a) 

-fiAqo + (Vixo)^9o - {uq ■ V)9o + Vtto = 0 in i7, (5.55b) 

divuo = 0, divQo = 0 in ill, (5.55c) 

^0 = 9, 9 o = 0 on tin n (5.55d) 

riQ = 9 o = 0 in B, (5.55e) 

rtg = 0, 90 - a, on T, (5.55f) 
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together with the following gradient equation: 


K 7 + 2Ao + fJ-di/Qo ■ di^UQ + / • a = 0 on F, 


(5.56) 


which is consistent with the adjoint system (2.7) and the strong form of (2.6) from [3], 
taking into account the volume constraint (see ( 4.32[ )) and the additional perimeter regu¬ 
larization. 


Remark 5.1 (Linear scaling for the correction constant 5^). Suppose 5^ = £, then we 


observe from (5.3) that 


i.e., 


2VV'(^o) 


Mo($o) = Mi(4>o)$i= ^ 

v2co v2co 


The presence of this extra factor of 2Wi(<f>o)‘l’i alters the jump term of (5.51) 

to 

[d,Qi ■ d,Ui - d,Ui ■ a(Alo(^o)9.^i + Mi(4>o)^i^o)]!: 


] Z=+oo 
Z=-oo 


a 


= [di,qo ■ di,uo]l - — 

2co 


^’n 




dzUi 


a 


— duQo ' di/iiQ ■ dijiiQ,. 

2co 


where we have used (5.29). Thus, instead of (5.56), we obtain 


Kj + 2Ao + pdi,qo ■ di,UQ + -^d^UQ ■ a + f ■ a = 0 on F. 

2co 


6 Numerical computations 


In this section we investigate the phase field approach numerically. We minimize the drag 
and maximize the lift-to-drag ratio of an obstacle in outer flow and apply both phase field 
approximations of the corresponding surface functionals. 

Concerning numerical results in the literature we refer to the minimization of the drag 
functional in [261 16]. where a sharp interface approach is used. In m the porous medium 
approach is used, where the authors argue, that the term asU^ is a valid approximation 
for the hydrodynamic force. 

Let us start with defining the free energy f:. Here we use 

Hv) = ^ (max^(0,y- 1) +min2(0,y+ 1)) -h J(1 

1 ( 6 . 1 ) 

^2(s-l)- 


V'(y) 




Note that can be obtained by using a Moreau-Yosida relaxation of the double¬ 
obstacle free energy (3.12) with the relaxation (or penalization) parameter s » 1, and the 
scaling of the argument and the shifting are chosen such that ip has its minima at y = ±1 
with "(/^(il) = 0. 
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We further introduce the convex-concave splitting 


Tp = 'lp++'lp_, 

V'+(2/) = ^ (max^ ^0, - 1 j 


+ min 


1 




2(s-l)’ 


where 'p+ is the convex part of and i/’- is its concave part. 
Next we define the interpolation function as 


ae{y) = “ 


0 


if ?/> 1 , 
if 1 > (f > 9, 


( 6 . 2 ) 


(i-e)(3+e)(y 

where a is a given constant, and we choose 6 = 0.99. This function a^{y) describes a 
linear function between y - -2 and y = 9 and has a quadratic extension between y = 9 
and y - 1. We fulhll Assumption 3.2 with Sa = -2 and s\, - 1. Note that we do not fulhll 


the regularity e (M) at Sa- But this is not a severe violation since in practice it 
holds that -2 < (p^ and we can control the violation of the bound -1 < ps by choosing an 
appropriate relaxation parameter s. 

For solving the optimization problem (3.8) we use a mass conserving Ff'^-gradient flow 


approach, following m- For this purpose we introduce an artificial time variable t and 
solve the following evolution equation for the phase field variable Pe{t) which is obtained 
from ( 4.30| ): 

dtPe = AWe, 


We = -le^Pe + '^'tp'iy^e) + a^i^e) ^ 


1 , 
2 ' 


'^e ' Qs\ Jip: 


(6.3) 


J^f, = - div {M{pe)f^4h{x,VUs,Pe,Vp£)) , 


where is obtained from (3.4), is obtained from (4.25) and abbreviates the terms 


arising from the differentiation of the functional h, as shown in Theorem 4.10, Note that 


we include the factor ^ into the parameter 7 . Using the gradient flow approach allows 
us to use nonlinear parts of the gradient, for example the derivative of implicitly in 
time in a time stepping scheme, which for the chosen free energy is favorable in view of 
stability reasons. 

After time discretization with variable time step size we at each time instance 
solve the following problem: 

Given find u^, p^, q^, and vTg fulfilling the primal system 


ae{p^)u£ - pAu£ + {u£ ■V)u£ + Vp£ = f in n, 

divUg = 0 in n, 

U£ = g on dfl, 

the adjoint system 

a£{p’p)qe - /idiv (VQe + (Vq^)^) + {VUef qe - (we ■ V)qe + VvTe 

= ae{pl)u£ - div [M{pl)T> 2 h) in 0 , 

divq^ =-M.{p^)Y)^h + §£ in fl, 

Qe = 0 on (9U, 


(6.4) 


(6.5) 
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and the Cahn-Hilliard system 


fc+1 Ai+1 A fc+1 , k 

'e = '^ + (p^ 

in H, 

y = -7eA9,**‘ + 1 



in H, 

X = M'{p^)h{x,VUe,Pe,Vp^'*'^) 


- div [M{p'^)I) 4 h{x, \/Ue,Pe, VPe^^)) , 


0 = ■ V'gfi + M{p’^)ugn ■ DX 

on 

0 = ■ vgn 

on d^l 


As noted above, we evaluate at the new time instance for stability reasons. 

For the spatial discretization piecewise linear and globally continuous finite elements 
are used for the variables Pe, and vTe, while piecewise quadratic and globally 

continuous elements are used for and q^- The meshes are adapted using the jumps of 
the normal derivative of and over edges of the underlying discretization mesh, 
see [a ED], together with a Dorfler marking |10j . 


6.1 Minimization of the hydrodynamic force of an obstacle 

We investigate the minimization of the drag of an obstacle of fixed area in a channel flow 
with block inflow profile. 

The computational domain is = (0,1.7) x (0,0.4). The initial phase field is defined 
as a circle of radius r = 0.05 with center at M = (0.5,0.2). The boundary velocity is set to 
g{x,y) = (1,0)^. We fix = 0, s = 1 x 10®, and / = 0. We further set 

:= ?nun(/iT||Vu;^||22(T)), 


where the minimization is carried out over all triangles T. Here, the diameter of triangle 
T is denoted by hx, and ^ is a positive scaling parameter typically set to ^ = 5. This 
CFL-like condition prevents the interfacial region from moving too fast for the adaptation 
process. 

We restate the definition of the phase field approximation of the hydrodynamic force 
in a direction a as 

:= J^M{ips)S/(ps ■ (l* (VWe + (VUg)^) - Pel) • adx. (6.7) 


When a is equal to the direction of the flow, i.e., a = (1,0)^, we denote the resulting 
approximation as , which corresponds to the drag of the obstacle. Meanwhile, if a is 
perpendicular to the direction of the flow, i.e., a = (0,1)^, then we denote the resulting 
approximation as F^, which corresponds to the lift of the obstacle. 

From (4.40) and (4.43), the terms arising from the derivatives of h in systems (6.5) 


and 


in the present setting are given as 


(-div [M{ip^)T>2h) ,v) = p ■ (Vu + (Vu)^) adx Vu e Hq{Q), 

(-A4(<p^)D3/i + 7?e,r?) = £^M{(p^)Vp^-adx^ridx \fr]€Ll{n), 

“ (^e ' V)^^) ■ < dx VC e H^(H). 
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Figure 1: Result for minimizing the drag using A4((pe) = the left plot we 

show the obstacle (i.e., (pe < 0 ) and streamlines of in black, and the pressure outside 
of the obstacle in gray. Darker gray means higher pressure. On the right we show lugl 
in gray, where darker gray means lower velocity. The isoline e 0 is shown in white 
and again streamlines are displayed in black. The results for A4((f^) = | are visually 
indistinguishable from these results. Note that we only show the computational domain 
in the neighbourhood of the obstacle. 


Next, we report on the numerical results for the case of minimizing F^. The param¬ 
eters are chosen as e = 0.00025, a = 0.03, p = 0.001, and 7 = 0.01. We note that we use 
path-following with respect to the value of /i, starting from fi = 0 . 01 , and also for the value 
of 7 , starting from 7 = 0.1. In Figure!^ we show results obtained with our approach. 

The drag for Ai((ps) = is given by = 3.9454 x 10'^ (3.9492 x 10'^), and 

for A4((pe) = I '''''6 have F^ - 3.9117 x 10“^ (3.9499 x 10“^). In brackets we give the drag 
obtained by evaluating the surface formulation over the isoline = 0. We see that both 
formulations give very similar results. 

6.2 Maximization of the lift-drag ratio of an obstacle 

Based on the results of the previous section we now investigate the maximization of the 
lift-to-drag ratio given by 

R:= F^jF^, 


To this end, we consider 



X, Vii£,P£, V(^£)dx 


+ (V^s)^) -p£l)a-^dx 
+ (VUs)'^) -Pel)adx ’ 


with a = ( 1 , 0 )^ and = ( 0 , 1 )^. 

The numerical setup is the same as in the previous section and the parameters are 
chosen as e = 0.0005, a = 4, p = 1/15, and 7 = 0.3. In this example we fix the y-coordinate 
of the center of mass of the obstacle by a Lagrange multiplier approach in order to keep 
it fixed at the initial position. We define the center of mass of the obstacle as 


com = 


/o^dx 


In Figurewe show results for this parameter set. 

We observe the expected optimal shape for both formulations, but for Ai{pe) = 

obtain a longer and thinner obstacle. 

The lift-to-drag ratio for Ai{pe) = is R - 1.1104, and for M.{pe) = | it is 

R = 0.9885. We stress that, here we calculate with a rather small value of p, = 1/15 and 
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Figure 2: Result for maximizing the lift-to-drag ratio using (left) and 

A4((fc) = I (right). The obstacle (i.e., < 0) and streamlines are shown in black and the 

velocity magnitude in gray. Darker gray means larger velocity. Note that we only show 
the computational domain in the neighbourhood of the obstacle. 


that the minimal magnitude of velocity inside the obstacle is 4x10 which is rather large. 
However, we think that the results are a promising starting point for further investigations. 


References 

[1] L. Ambrosio, N. Fusco, and D. Pallara. Functions of Bounded Variation and Free 
Discontinuity Problems (Oxford Mathematical Monographs). Oxford University Press, 
USA, 2000. 

[2] J. A. Bello, E. Fernandex-Cara, J. Lemoine, and J. Simon. The differentiability of 
the drag with respect to the variations of a Lipschitz domain in a Navier-Stokes flow. 
SIAM Journal on Control and Optimization, 35(2):626-640, 1997. 

[3] S. Boisgerault and J.P. Zolesio. Shape derivative of sharp functionals governed by 
Navier-Stokes flow. In W. Jager, J. Necas, O. John, K. Najzar, and J. Stara, ed¬ 
itors, Partial Differential Equations: Theory and Numerical Solution, pages 49-63. 
Chapman and Hall/CRC, 1993. 

[4] T. Borrvall and J. Petersson. Topology optimization of fluids in Stokes flow. Internat. 
J. Numer. Methods Fluids, 41(1):77-107, 2003. 

[5] B. Bourdin and A. Chambolle. Design-dependent loads in topology optimization. 
ESAIM Control Optim. Calc. Var., 9:19-48, 8 2003. 

[6] C. Brandenburg, F. Lindemann, M. Ulbrich, and S. Ulbrich. A Continuous Adjoint 
Approach to Shape Optimization for Navier Stokes Flow. In K. Kunisch, J. Sprekels, 
G. Leugering, and F. Troltzsch, editors. Optimal Control of Coupled Systems of Par¬ 
tial Differential Equations, volume 158 of Internat. Ser. Numer. Math., pages 35-56. 
Birkhauser, 2009. 

[7] C. Brandenburg, F. Lindemann, M. Ulbrich, and S. Ulbrich. Advanced Numerical 
Methods for PDE Constrained Optimization with Application to Optimal Design in 
Navier Stokes Flow. In G. Leugering, S. Engell, A. Griewank, M. Hinze, R. Rannacher, 
V. Schulz, M. Ulbrich, and S. Ulbrich, editors, Constrained Optimization and Optimal 
Control for Partial Differential Equations, pages 257-275. Birkhauser, 2012. 

[8] C. Carstensen and R. Verfiirth. Edge Residuals Dominate A Posteriori Error Esti¬ 
mates for Low Order Finite Element Methods. SIAM J. Numer. Anal., 36(5): 1571- 
1587, 1999. 


38 





































[9] X. Chen. Global asymptotic limit of solutions of the Cahn-Hilliard equation. Journal 
of Differential Geometry, 44(2):262-311, 1996. 

[10] W. Dorfler. A convergent adaptive algorithm for Poisson’s equation. SIAM J. Numer. 
Anal, 33(3):1106-1124, 1996. 

[11] L.C. Evans and R.F. Gariepy. Measure Theory and Fine Properties of Functions. 
Studies in advanced mathematics. GRC Press, 1992. 

[12] P.G. Fife and O. Penrose. Interfacial dynamics for thermodynamically consistent 
phase-field models with nonconserved order parameter. EJDE, 16:1-49, 1995. 

[13] G.P. Galdi. An Introduction to the Mathematical Theory of the Navier-Stokes Equa¬ 
tions. Springer Monographs in Mathematics. Springer, New York, NY, 2011. 

[14] H. Garcke and G. Hecht. Applying a phase field approach for shape optimization of a 
stationary Navier-Stokes flow, to appear in ESAIM: Control Optim. Calc. Var., 2014. 

[15] H. Garcke, G. Hecht, M. Hinze, and G. Kahle. Numerical approximation of phase 
field based shape and topology optimization for fluids. arXiv: 1405.3480, 2014. 

[16] H. Garcke and B. Stinner. Second order phase field asymptotics for multi-component 
systems. Interfaces and Free Boundaries, 8:131-157, 2006. 

[17] E. Giusti. Minimal surfaces and functions of bounded variation. Notes on pure 
mathematics. Dept, of Pure Mathematics, 1977. 

[18] C. Hecht. Shape and topology optimization in fluids using a phase field approach and 
an application in structural optimization. Dissertation, University of Regensburg, 
2014. 

[19] M. Hinze, R. Pinnau, M. Ulbrich, and S. Ulbrich. Optimization with PDE constraints. 
Mathematical Modelling: Theory and Applications, Volume 28. Springer, 2008. 

[20] B. Kawohl, A. Cellina, and A. Ornelas. Optimal Shape Design: Lectures Given at 
the Joint C.I.M./C.I.M.E. Summer School Held in Troia (Portugal), June 1-6, 1998. 
Lecture Notes in Mathematics / C.I.M.E. Foundation Subseries. Springer, 2000. 

[21] T. Kondoh, T. Matsumori, and A. Kawamoto. Drag minimization and lift maximiza¬ 
tion in laminar flows via topology optimization employing simple objective function 
expressions based on body force integration. Structural and Multidisciplinary Opti¬ 
mization, 45(5):693-701, 2012. 

[22] L. Modica. The gradient theory of phase transitions and the minimal interface crite¬ 
rion. Arch. Ration. Mech. Anal, 98(2):123-142, 1987. 

[23] F. Murat. Gontre-exemples pour divers problemes ou le controle intervient dans les 
coefficients. Annali di Matematica Pura ed Applicata, 112(1):49-68, 1977. 

[24] O. Pironneau. On optimum design in fluid mechanics. J. Fluid Mech., 64:97-110, 5 
1974. 

[25] P.I. Plotnikov and J. Sokolowski. Shape derivative of drag functional. SIAM J. 
Control Optim., 48(7):4680-4706, 2010. 


39 



[26] S. Schmidt and V. Schulz. Shape Derivatives for General Objective Functions and the 
Incompressible Navier-Stokes Equations. Control Cybernet., 39(3):677-713, 2010. 

[27] H. Sohr. The Navier-Stokes Equations: An Elementary Functional Analytic Approach. 
Birkhauser Advanced Texts. Springer Verlag, 2001. 

[28] L. Tartar. Problemes de Controle des Coefficients Dans des Equations aux Derivees 
Partielles. In A. Bensoussan and J.L. Lions, editors, Control Theory, Numerical 
Methods and Computer Systems Modelling, volume 107 of Lecture Notes in Economics 
and Mathematical Systems, pages 420-426. Springer, 1975. 

[29] F. Troltzsch. Optimal Control of Partial Differential Equations: Theory, Methods, 
and Applications. Graduate studies in mathematics. American Mathematical Society, 
2010. 

[30] R. Verfiirth. A review of a posteriori error estimation and adaptive mesh-refinement 
techniques. Wiley-Teubner series; Advances in Numerical Mathematics. Wiley- 
Teubner, New York, 1996. 

[31] E. Zeidler. Nonlinear Functional Analysis and Its Applications: Part I: Fixed-Point 
Theorems. Springer, New York, Berlin, 1986. 

[32] E. Zeidler. Nonlinear Functional Analysis and Its Applications: Part IV: Applications 
to Mathematical Physics. Springer Verlag, New York, 1997. 


40 



